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PhEFACE 


A technical  workshop  on  viscous  interaction  and  boundary- 
layer  separation  was  held  at  The  Ohio  State  University  on  August 
lo-17,  197o.  The  purpose  of  the  workshop  was  to  assess  the  cur- 
rent state  of  progress  and  to  stimulate  further  efforts  toward 
the  Foal  of  practical  and  reliaule  analysis  of  viscous  flows  at 
high  Reynolds  numbers.  For  definiteness , the  scope  of  the  work- 
shop was  limited  to  theoretical  aspects  of  the  subject,  with  em- 
phasis placed  on  fundamental  ideas  and  problem  areas  rather  than 
on  details  of  computations  and  results. 

Only  abstracts  of  the  papers  were  requested  in  order  to  pro- 
mote informality  and  a free  interchange  of  ideas.  Each  participant 
was  given  twenty  minutes  for  a presentation  on  recent  research 
efforts  (in  most  cases  work  in  progress  was  reported),  followed 
oy  approximately  ten  minutes  of  lively  discussion.  The  workshop 
was  opened  by  the  keynote  lecture  presented  by  Professor  Keith 
Stewartson  on  'Sychev's  Theory  of  Separation'.  In  addition,  an 
evening  session  was  devoted  to  unstructured  discussion. 

Acknowledgement  is  due  Capt . Scott  McKae , AFFDL,  and  Lt . Col. 
Robert  Smith,  AFOSR,  for  their  assistance  in  obtaining  the  funds 
that  made  the  workshop  possible. 
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INTRODUCTION 


The  collected  abstracts  of  the  workshop  lectures  are  presented 
in  these  Proceedings.  In  addition,  the  Organizing  Committee  thought 
it  worthwhile  to  delineate  some  of  the  more  significant  points  of 
the  lectures  and  discussions.  These  fall  naturally  into  several 
more  or  less  distinct  groupings: 

Numerical  problems; 

Separation  bubbles; 

Breakaway  separation; 

Trailing  edges  (steady  and  unsteady); 

Tubes  and  channels; 

Three-dimensional  interactions. 

A recapitulation  of  the  main  points  brought  out  in  the  discussion 
follows  under  each  of  the  above  group  headings.  It  may  be  remarked 
that  some  of  this  material  is  controversial  even  within  the  Organ- 
izing Committee. 


Numerical  Problems 

Approximate  analytical  theories  are  not  as  useful  to  the  engi- 
neer as  formerly  because  of  the  near  universal  availability  of  large 
scale  computers  and  attendant  software  for  solving  problems  in  fluid 
dynamics.  Nevertheless,  numerical  methods  are  not  without  problems, 
and  anlytical  solutions  are  useful  not  only  in  understanding  the 
flow  processes  but  also  as  a guide  to  resolving  the  numerical  diffi- 
culties. 

Tne  numerical  problems  which  one  encounters  in  high  Reynolds 
number  separated  flows  fall  into  several  categories.  First,  the 
asymptotic  analysis  has  revealed  that  at  high  Reynolds  number, 
severe  scaling  problems  exist  at  and  downstream  of  separation. 

Tne  scalings  indicated  by  the  asymptotic  analysis  have  generally 
not  been  honored  in  Navier-Stokes  solution  methods  to  date  and 
tnerefore  there  is  some  doubt  about  the  accuracy  of  many  such  sol- 
utions and  the  reliability  of  the  method.  None-the-less , some 
Navier-Stokes  solutions  have  shown  good  agreement  with  available 
experimental  data,  even  though  no  direct  effort  was  made  to  honor 
the  asymptotic  scale  laws.  These  scalings  have  been  honored  in 
some  boundary- layer-like  calculations  and  the  results  of  these  cal- 
culations tend  to  confirm  the  correctness  of  the  triple-deck  analy- 
sis . 

Second,  now  that  the  nature  of  interaction  and  separation  is 
more  clearly  understood,  it  would  appear  that  this  knowledge  would 
be  useful  in  formulating  numerical  techniques.  This  has  been  done 
in  boundary-layer  solution  methods,  but  little  has  been  done  in 
this  direction  in  terms  of  solution  of  the  Navier-Stokes  equations. 

I 


1 

I 


Third,  it  would  appear  that  significant  improvements  can  still 
be  made  in  the  computational  convergence  rate  of  numerical  methods 
for  separated  flows  even  when  the  flow  is  governed  by  the  simplest 
equations  possible,  the  lower  deck  equations.  In  the  subsonic 
case,  for  example,  the  convergence  of  present  methods  is  extremely 
slow  and  requires  severe  under-relaxation.  Some  mixed  inverse- 
direct  iteration  methods  show  promise  in  this  area.  At  the  same 
time  it  would  appear  that  work  still  needs  to  be  done  on  developing 
better  difference  and  switching  methods  for  reverse-flow  regions. 

Finally,  we  still  need  to  resolve  many  remaining  questions 
about  whether  implicit  or  explicit,  conservative  or  non-conservative, 
splitting  or  non-splitting,  methods  are  most  efficient  for  solving 
the  Navier-Stokes  equations  at  high  Reynolds  number.  For  example, 
some  believe  that  a problem  can  be  solved  most  efficiently  with 
totally  implicit  methods  while  others  believe  that  the  viscous 
portion  of  the  flow  should  be  done  with  an  implicit  method  and 
the  inviscid  portion  done  with  an  explicit  method. 


Closed  Separation  Bubbles 

It  appears  that  a reliable  theoretical  model  is  now  available 
for  the  case  of  small  two-dimensional  laminar  separation  bubbles 
in  steady  flow.  This  model  is  based  on  the  new  asymptotic  limit 
solution  obtained  using  a three-deck  model  of  the  interaction/ 
separation  process.  It  remains  now  to  seek  ways  of  solving  the 
resulting  governing  equations  as  efficiently  as  possible. 

The  triple-deck  model  provides  a sound  theoretical  basis  for 
the  use  of  the  interacting  boundary-layer  equations  which  can  be 
viewed  as  a single  set  of  governing  equations  uniformly  valid 
across  the  bottom  two  decks  of  the  triple-deck  model.  Evidence 
to  support  the  use  of  the  interacting  boundary-layer  approach  for 
steady  laminar  flows  continues  to  amass  for  a wide  variety  of  flow 
configurations  and  flight  regimes.  The  extension  of  this  approach 
to  the  turbulent  flow  case  has  shown  an  equal  level  of  success  al- 
though, as  yet,  a helpful  asymptotic  model  has  not  been  developed. 

The  principle  difficulties  remaining  for  closed  separation 
bubbles  are  those  concerned  with  either  developing  better  methods 
of  solving  the  currently  available  model  equations,  or  extending 
the  present  models  to  increase  their  utility.  The  following  is  a 
list  of  some  of  the  problems  in  the  first  category:  (a)  solution 

of  the  highly  nonlinear  displacement-interaction  equations  pre- 
sently requires  use  of  inefficient  local  interaction  techniques; 

(b)  the  global  interaction  effect  (downstream  to  upstream)  of  sub- 
sonic flows  currently  involves  slow  iteration  techniques;  (c)  large 
separation  bubbles  develop  an  inviscid  core  for  which  current  me- 
thods are  inefficient;  (d)  the  penetration  of  shock  waves  deep 
into  separating  supersonic  turbulent  and  hypersonic  laminar  bound- 
ary layers  obliterates  numerical  accuracy  in  current  approaches; 
and  (e)  the  scale  of  separation  processes  are  very  small  and  diffi- 
cult to  resolve  with  current  techniques. 
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For  the  second  category  of  problems,  those  involving  analy- 
tical extension  of  current  models,  there  are  currently  three  pro- 
blems that  need  attention;  (a)  the  inclusion  of  imbedded  shock 
waves  in  the  supersonic/hypersonic  separation  model;  (b)  the  ex- 
tension of  asymptotic  concepts  to  the  turbulent  separation  problem; 
and  (c)  the  development  of  an  extended  limit  solution  for  the 
supersonic  separation  case  to  allow  application  over  a more  prac- 
tical range  of  parameters. 


Breakaway  Separation 

Bluff  bodies  and  airfoils  at  high  angle  of  attack  are  charac- 
terized by  breakaway  separation  for  which  the  flow  does  not  reattach 
on  the  body.  The  problem  is  an  old  one  in  the  literature  of  fluid 
mechanics,  and  the  paradoxes  associated  with  the  classical  Helmholz- 
Kircnhoff  free-streamline  theory  are  well  known.*  Briefly,  the  in- 
viscid  theory  requires  the  separation  point  to  occur  in  a region  of 
favorable  pressure  gradient  on  a smooth  bluff  body,  contradicting 
the  ideas  of  classical  boundary-layer  theory.  The  paradox  appears 
to  be  resolved  by  Sychev's  triple-deck  theory  involving  interaction 
between  the  viscous  and  inviscid  fluid  layers;  in  this  model  the 
triple-deck  region  with  adverse  pressure  gradient  for  separation  is 
embedded  in  the  potential  flow  region  of  favorable  pressure  gra- 
dient. Up  to  the  time  of  the  workshop,  it  was  uncertain  as  to 
whether  Sychev's  theory  is  well-founded  or  not,  since  no  solution 
of  his  equations  had  yet  been  obtained.  However,  one  speaker  did 
report  on  progress  toward  a solution  with  results  suggesting  that 
a solution  exists  and  is  unique  (F.T.  Smith;  work  in  progress,  not 
described  in  abstract). 

It  is  clear  that  more  theoretical  work  is  needed  along  these 
lines;  for  example,  treatment  of  the  case  of  the  turbulent  boundary 
layer,  and  inclusion  of  unsteady  effects  that  occur  at  large  Rey- 
nolds number  even  in  the  laminar  case.  Nevertheless,  certain  pro- 
perties of  the  present  theory  may  be  useful  in  guiding  numerical 
analyses  of  bluff-body  flows.  In  particular,  the  pressure  distri- 
bution involves  a rather  rapid  rise  In  the  immediate  vicinity  of 
the  separation  point,  requiring  a rather  fine  mesh  to  reproduce 
by  finite  differences,  together  with  a slow  decay  over  a large 
streamwise  distance  to  the  "outer"  inviscid  flow  both  upstream 
and  downstream,  thus  requiring  the  numerical  mesh  to  extend  over 
a large  region  about  the  separation  point.  These  opposing  trends 
clarify  the  well-known  difficulty  of  computing  accurate  solutions 
of  bluff-body  flows  at  large  Reynolds  number. 

As  a final  point,  there  is  suggested  the  possibility  of  in- 
corporating these  triple-deck  concepts  into  ad  hoc  analyses  of, 
for  example,  separated  wakes  of  airfoils  at  high  angle  of  attack. 


* See  for  example,  G.  Birkhoff,  Hydrodynamics , a Study  in  Logic , 
Fact , and  Similitude , Princeton  University  Tress  ( 1950 ) . 
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Trailing  Edges 

Viscous  interaction  at  a trailing  edge  results  from  the  abrupt 
change-over  from  the  no-slip  condition  on  the  airfoil  surface  to 
the  wake-symmetry  condition  downstream.  For  the  flat  plate  in 
incompressible  flow  at  zero  incidence,  the  pressure  falls  and  the 
sKin  friction  rises  (to  a value  nearly  35$  greater  than  the  Blasius 
value)  as  the  flow  is  accelerated  toward  the  trailing  edge.  Down- 
stream the  pressure  rises  very  rapidly,  overshoots  the  freestream 
value,  and  then  decays  very  slowly  over  a large  distance.  These 
results  from  the  triple-deck  theory  have  strong  implications  when 
choosing  mesh  distributions  for  solving  the  Navier-Stokes  equa- 
tions numerically.  An  even  finer  substructure  than  the  triple  deck 
was  described  at  tne  workshop;  however,  this  flow  structure  is  pas- 
sively driven  by  that  in  the  triple-deck  and  consequently  does  not 
alter  tne  larger  scale  properties. 

The  proulem  of  trailing-edge  flow  on  an  airfoil  with  angle-of- 
attack  has  now  been  treated  in  the  framework  of  triple-deck  theory, 
both  for  subsonic  and  supersonic  freestream.  The  most  obvious 
practical  application  is  the  reduction  of  lift  due  to  viscous  in- 
teraction through  weakening  of  tne  Kutta  condition.  The  flow  has 
been  computed  for  angle-of-attack  up  to  the  point  of  traii ing-edge 
separation.  For  larger  angles  the  flow  structure  is  not  known, 
obviously  more  work  is  needed  on  this  problem. 

For  turbulent  flow,  a three-deck  asymptotic  flow  structure 
exists  in  the  limit  of  infinite  Reynolds  number,  somewhat  like  the 
laminar  triple-deck,  but  different  in  that  the  streamwise  extent 
is  short  compared  with  the  boundary-layer  thickness.  Experiments 
were  reported  in  which  separation  bubbles  were  observed  as  short 
as  20$  of  the  boundary-layer  tnickness.  This  short  length  scale, 
over  which  large  pressure  changes  occur,  again  imposes  severe 
restrictions  on  a turbulent  "Navier-Stokes"  computation.  As  for 
tne  laminar  case,  the  flow  structure  beyond  the  stall  point  is  not 
yet  known. 

The  unsteady  flow  structure  was  outlined  for  oscillatory  mo- 
tion. The  dominant  effect  is  the  unsteady  triple-deck  at  the  trail- 
ing edge,  which  produces  unsteady  contributions  to  the  lift  and 
pitching  moment  having  phase  shifts  of  A 5 degrees  and  90  degrees. 
Such  phase  shifts  could  have  a significant  influence  on  flutter 
boundaries.  Much  more  work  is  needed  in  this  area;  for  example, 
numerical  solutions  of  tne  unsteady  triple-deck  equations  have  not 
yet  been  obtained. 


and  Channel: 


Numerical  solutions  to  parabolized  tube  and  channel  flow  equa- 
tions have  indicated  the  possibility  of  branching  behavior  in  the 
numerical  solutions  as  the  streamwise  step  size  is  decreased.  Tills 
is  analagous  to  the  behavior  which  lias  been  previously  found  in 
supersonic  interacting,  boundary  layers  on  flat  plates  for  example. 


It  has  been  found  that  as  in  the  case  of  an  external  flow,  a triple- 
dec*c  analysis  is  applicable  to  channel  and  tube  flows,  but  the  ex- 
pansions are  different.  The  asymptotic  analysis  reveals  an  upstream 
influence  and  therefore  provides  a mechanism  for  the  branching  be- 
havior observed  numerically.  More  work  needs  to  be  done  to  demon- 
strate the  connection  between  the  asymptotic  theory  and  the  observed 
urancning  behavior  of  the  numerical  solutions. 

For  curved  tubes  the  far  downstream  solution  is  not  known.  A 
possible  formulation  of  the  appropriate  asymptotic  model  and  a status 
report  on  the  numerical  solution  of  the  resulting  equations  was 
given . 


Three-dimensional  Interactions 


Almost  all  work  in  viscous-interaction  has  been  carrried  out 
on  two-dimensional  flow  problems  (we  include  axi-symmetric  flows 
in  this  category).  Indeed,  all  but  two  of  the  abstracts  presented 
for  this  workshop  were  concerned  with  two-dimensional  flows.  A 
three-dimensional  version  of  the  triple-deck  theory  was  described 
and,  as  one  application,  similarity  laws  were  derived  for  pressure 
and  shear-stress  distributions  on  simple  swept  configurations. 
Preliminary  results  were  presented  as  well  for  more  general  three- 
dimensional  interacting  flows.  Clearly  this  field  contains  a wide 
variety  of  problems  for  future  applications. 


PROGRAM 

WORKSHOP  ON  VISCOUS  INTERACTION  AND 
BOUNDARY-LAYER  SEPARATION 


Monday  horning  (August  16) 

Session  1 (8:30  A.M. ) 

K.  Stewartson  - Sychev  Theory  of  Separation. 

S.  WeinDaum  - A New  Pressure  Hypothesis  and  Approximate 

Theory  for  Streaming  Motion  Past  Bluff  Bodies  at 
Re  from  0(1)  to  0(102). 

J.  Carter  - An  Inverse  Technique  for  Viscous-inviscid 
Interaction . 

Coffee  Break  (10:05  A.M.) 

Session  2 (10:20  A.M.) 

R.  MacCormack  - An  Efficient  Numerical  Method  for  Solving 
Time-dependent  Compressible  Navier-Stokes  Equations 
at  High  Reynolds  Number. 

J.  Rom  - The  Calculation  Methods  for  Flows  with  Strong 

Interactions  Between  Viscous  and  Inviscid  Regions. 

T.  Adamson  - Shock-wave-turbulent  Boundary  Layer  Inter- 

action at  Incipient  Separation. 

A.  Nliiwick  - Freely  Interacting  Transonic  Boundary  Layers. 

Lunch  (12:30  P.M. ) 

Monday  Afternoon  (August  16) 

Session  1 (2:00  P.M.) 

N.  Riley  - Flow  Separation  from  a Smooth  Surface  in  a 
Slender  Conical  Flow. 

F.  Blottner  - Viscous  Slender  Channel  Flows. 

J.-C.  Le  Balleur  - Viscous  Interaction  Including  Turbulent 
Boundary  Layer  Separation  in  Two-Dimensional  Flow. 

P.  Roache  - Difficulties  with  Brute-force  Numerical  Solu- 
tions of  Navier-Stokes  Equations  for  High  - Re  Separ- 
ated Flows . 

Coffee  Break  (3:J5  P.M.) 

Session  2 (3:50  P.M.) 

R.  Melnik  - Turbulent  Interaction  at  Trailing  Edges. 

S.  Brown  - Free  Interactions  in  Unsteady  Boundary  Layers 

M.  Werle  - numerical  Studies  of  High  Reynolds  Number  Laminar 
Separated  Flow. 

G.  Kuhn  - Calculation  of  Axisymmetric  Separated  Turbulent 

Boundary  Layers. 

Dinner  (7:00  P.M.) 


b 


I 

I 


1 


Tuesday  Morning  (August  17) 


Session 

S. 

F. 
W . 

A. 


1 (8:30  A.M.  ) 

Dennis  - Boundary-layer  Calculations  for  a Circular 
Cylinder. 

Smith  - Tube-flows:  dn  the  effects  of  Constrictions, 

Branchings , upstream  Influence  and  Separation. 

Hankey  - A Controlling  Factor  in  the  establishment 
of  the  Size  of  tne  reparation  Region. 

van  de  Vooren  - Viscous  Flow  Hear  the  Trailing  Edge 
of  a Flat  Plate. 


Coffee  Break  (10:05  A.M. ) 


Session  2 (10:20  A.M.) 

rl.  K.  Cheng  - Study  on  Unsteady  Trai ling-edge  Stall. 

H.  McDonald  - Viscous  Interaction  and  Separation. 

P.  Williams  - Reversed  Flow  in  Boundary  Layers  with  Prescribed 
Displacement  Thickness. 

R.  Ackerberg  - Separation  at  a Free  Streamline. 

Luncn  (12:30  P . M. ) 

Tuesday  Afternoon  (August  17) 


Session 

T. 

M. 

G. 

O • 

0. 

s. 


1 (2:00  P.M. ) 

Cebeci  - Separated  Flows  and  Their  Representation  by 
Boundary -Layer  Equations. 

Holden  - Experimental  Investigation  of  Current  Problems 
Involving  Boundary  Layer  Separation. 

Inger  - Upstream  Influence  in  Interacting  Mon-separated 
Turbulent  Boundary  Layers. 

Dennis  and  N.  Riley  - Flow  in  a Curved  Pipe  at  High  Dean 
Number. 

Burggraf  - The  Three-dimensional  Triple  Deck. 

Rubin  - Viscous  Interactions  in  Three  Dimensions,  (paper 
not  presented  due  to  illness.) 
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SEPARATION  AT  A FREE  STREAMLINE 


R.C.  Ackerberg 

Polytechnic  Institute  of  New  York 
Department  of  Chemical  Engineering 
Brooklyn,  New  York  11201 

ABSTRACT 

This  talk  will  deal  with  the  separation  that  occurs  in  the 
neighborhood  of  a sharp  trailing  edge  to  which  a free  streamline 
is  attached.  Previous  results  '2'  will  be  summarized,  and  some 
recent  work  dealing  with  large  capillary  effects  will  be  discussed. 

REFERENCES 

1.  Ackerberg,  R.C.,  " Boundary- layer  separation  at  a free  stream- 
line. Part  1.  Two-dimensional  flow,"  J.  Fluid  Mech.,  44, 
pp.  211-225,  1970. 

2.  Ackerberg,  R.C.,  "Boundary- layer  separation  at  a free  stream- 
line. Part  3.  Axisymmetric  flow  and  the  flow  downstream  of 
separation,"  J.  Fluid  Mech.,  59,  pp.  645-663,  1973. 

3.  Ackerberg,  R.C.,  "The  effects  of  capillarity  on  free  stream- 
line separation,"  J.  Fluid  Mech.,  _70,  pp.  333-352,  1975. 


SHOCK  WAVE-TURBULENT  BOUNDARY  LAYER 
INTERACTION  AT  INCIPIENT  SEPARATION 

T.  C.  Adamson,  Jr. 

The  University  of  Michigan 


ABSTRACT 


Asymptotic  methods  are  used  to  describe  the  structure  of  the  interaction 
region  formed  when  a shock  wave,  in  transonic  flow,  impinges  upon  a turbulent 
flat  plate  boundary  layer.  The  case  chosen  for  study  is  that  of  incipient  separa- 
tion. In  a previous  paper'  it  was  shown  that  just  as  in  the  weaker  shock  cases,2*3 
the  interaction  region  in  this  strong  shock  case  could  be  divided  into  regions 
such  that  the  flow  in  the  outer  region  (in  the  direction  normal  to  the  wall)  is 
governed  by  inviscid  flow  equations,  while  in  the  near  wall  regions  (Reynolds 
stress  sublayer  and  wall  layer)  Reynolds  and  viscous  stresses  must  be  taken 
into  account.  Unlike  the  weak  shock  cases,  there  are  also  two  regions  in  the 
flow  direction,  one  upstream  (inner)  region  with  characteristic  length  ordered 
by  the  extent  of  the  upstream  influence  and  thus  depending  upon  the  distance 
from  the  wall  to  the  undisturbed  sonic  line,  and  one  downstream  (outer)  region 
with  characteristic  length  depending  upon  the  boundary  layer  thickness.  Be- 
cause the  flowr  is  unseparated,  it  is  possible  to  find  solutions  in  the  inviscid  flow 
regions  independently  of  the  near  wall  regions,  and  these  solutions  lead  to  a 
relation  for  the  wall  pressure  distribution.  The  terms  in  the  w'all  pressure 
solution  are  then  used  in  the  wall  layer  solutions  to  find  the  relation  for  the  wall 
shear  stress. 

The  structure  of  the  near  w'all  regions  and  the  derivation  of  the  relation 
for  the  wall  shear  stress,  tw,  are  reviewed.  The  dependence  of  tw  upon  the 
closure  model  employed  is  discussed  as  is  the  variation  of  t with  distance 
along  the  wall  for  both  oblique  and  normal  impinging  shock  waves.  The  use  of 
the  expression  for  tw  in  forming  a separation  criterion  is  illustrated  and  the 
problems  associated  with  such  a calculation  are  discussed. 


REFERENCES 


1.  Adamson,  Jr.,  T.C.,  "Normal  Shock  Wave-Turbulent  B ounda  ry  La  ye  r Inte  r- 
actions  in  Transonic  Flow  Near  Separation,"  Workshop  on  Transonic  Prob- 
lems in  Turbomachinery,  Monterey,  California,  Feb.  1 1 -13,  1976.  Pro- 
ceedings to  be  published. 

L.  Adamson,  Jr.  , I . and  Feo,  A.  , "Interaction  Between  a Shock  Wave  and  a 
Turbulent  Boundary  Layer  in  Transonic  Flow,  " SIAM  J.  Appl.  Math.,  J9, 
pp.  Ill - 145,  1 975. 

3.  Melnik,  R.E.  and  Grossman,  B.  , "Analysis  of  the  Interaction  of  a Weak 
Normal  Shock  Wave  with  a Turbulent  Boundary  Layer,"  AIAA  Paper  No. 
74-598,  1974. 


VISCOUS  INTERACTION  INCLUDING  TURBULENT  BOUNDARY 
LAYER  SEPARATION  IN  TWO  DIMENSIONAL  FLOWS 


by  Jean-Claude  Le  Balleur 

Office  National  d' Etudes  ct  de  Recherches  Ac-rospntiales  (ONERA) 

92}20  CHATILLON  (France) 


A pc  tract  - 


A study  is  developped  to  provide  a numerical  prediction  of  steady  two- 
cimensional  flows  with  strong  viscous  interaction,  by  meansof  a mutched  viscid- 
inviscid  approximation. 

In  a conventional  way,  1’randtl  equations  are  used  to  describe  the 
turbulent  boundary  layer.  Recent  developments  on  laminar  separation  allow  to 
think  that  such  a simplification  is  consistent  with  incipient  separation  or 
small  separated  nubbles.  Because  of  practical  considerations,  eauations  arc 
s ->lvea  by  a moment  integral  method. 

As  usual  for  these  calculations,  the  matching  consists  in  searching 
common  boundary  conditions  at  the  outer  edge  of  the  dissipative  layer,  thus 
including  both  ordinary  boundary  layer  and  self-induced  separation.  Clip 
condition  on  displacement  surface,  or  injection  condition  at  the  wall  are  still 
equivalent,  and  the  more  convenient  is  used. 

Because  of  the  boundary  value  nature  of  the  matched  problem,  emphasis 
has  beer;  given  to  iterative  processes  where  repeated  inviscid  and  boundary 
layer  calculations  are  performed. 

A linear  analysis  first  explains  that  the  classical  iteration  on  the 
displacement  thickness  often  diverges,  but  however  provides  a suitable  method 
if  a locally  optimal  underrelaxation  is  added.  It  appears,  secondly,  for 
supersonic  regime,  that  the  method  correctly  includes  the  propagation  of  up- 
stream influence  as  long  as  the  boundary  layer  is  subcritical.  It  is  then 
nossible,  as  in  the  Werle  and  Vatsa  method,  to  prescribe  a downstream  pressure 
condition,  or  to  respect  automatically  a downstream  critical  point  condition, 
without  the  difficulties  associated  with  an  ill-posed  initial  value  problem. 
Thirdly,  convergence  of  the  matching  process  and  of  the  inviscid  algorithm 
may  be  simultaneously  performed  to  avoid  large  computing  times. 

This  method  failed  at  separation  because  of  the  boundary  layer  calculation. 
However,  following  Klineberg,  inverse  calculation  is  possible,  and  a linear 
analysis  shows  that  iteration  on  the  oressure  may  be  adequatly  underrelaxed 
to  ensure  convergence.  Despite  this,  inverse  boundary  layer  calculation  is  in- 
cot. venient  for  ordinary  attached  boundary  layer,  and  is  singular  at  the  critical 
pi int.;.  Consequently  inverse  matching  must  be  limited  to  separated  regions, 
a no  mixed  direct-inverse  matching  is  needed. 
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Complexity  and  practical  numerical  difficulties  led  us  tc  imagine  a 
new  method  where  only  the  boundary  layer  is  calculated  by  a direct  or  inverse 
method,  whereas  the  inviscid  calculation  is  always  a direct  one.  When  separa- 
tion occurs,  the  method  provides  an  automatic  change  of  the  displacement 
surface  with  regard  to  viscid  and  inviscid  estimations  of  the  pressure,  until 
convergence.  As  without  separation,  supersonic  upstream  influence  is  taken 
into  account  and  simultaneous  inviscid  resolution  - matching  iteration  is 
possible . 


The  method  seems  very  general  and  applies  form  subsonic  to  supersonic 
regime  to  any  inviscid  flow.  Qualitative  experimental  features  of  vi3cous 
effects  are  found  for  symmetric  arc-circular  transonic  airfoil  with  a small 
transonic  perturbation  inviscid  calculation,  and  for  supersonic  ramp  with 
Prandtl-Veyer  law. 


Koch 
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VISCOUS  SLlMMiR  C11ANNF.L  FLOWS 


F.G.  Blottner 
Sandia  Laboratories 
A1 buquerque , Mew  Mexico,  87115,  U.S.A. 

ABSTRACT 


ihe  slender  cbannel  approximation  lias  been  used  successfully  by  several 
authors  to  predict  the  flow  for  straight  nozzle  geometries  and  to  obtain  good 
agreement  with  experiments.  A finite-difference  scheme  and  computer  code  have 
been  developed!  for  solving  the  slender  channel  equations  in  curved  channels 
of  varying  height.  The  governing  equations  are  the  incompressible  boundary 
layer  equations  with  longitudinal  curvature  terns  included;  the  velocity  compo- 
nents at  the  wall  arc  zero  from  the  boundary  conditions.  New  independent  and 
dependent  variables  are  introduced  into  the  governing  equations  which  allow 
similar  solutions  to  be  obtained  and  the  computational  domain  becomes  a rect- 
angle. The  resulting  transformed  continuity  and  two  momentum  equations  are 
solved  with  an  implicit  second-order  difference  scheme  which  requires  the  solu- 
tion of  block  tridiagonal  difference  equations.  The  solution  procedure 
determines  cither  similar  solutions  which  can  be  used  as  initial  conditions  or 
nonsimilar  solutions  of  the  slender  channel  equations. 

Numerical  solution  of  the  flow  in  a converging  channel  have  been  obtained 
with  various  step-sizes  to  verify  the  accuracy  of  the  code.  The  flow  has  also 
been  calculated  in  a straight  expanding  channel  that  is  increasing  in  height  for 
a distance  S-r  and  then  is  of  constant  height  !L-  = (n^/n  )/-.  The  nond i mens i ona  1 
channel  height  variation  is  c c 

ne/n0  = t (1  + Hf)  + I (1  - Hf ) cos(Tir/4ST) 

where  = Ts/iTq,  s is  distance  along  the  channel  centerline  and  nG  is  initial 
channel  height.  When  S-]  = 10  and  Ilf  = 2,  there  is  a small  region  of  reverse 
flow.  Flie  separation  is  regular  and  no  numerical  problems  appear  to  occur. 
However  for  more  severe  separation  and  a small  step-size  in  the  marching  direc- 
tion, numerical  difficulties  arc  encountered  which  can  disappear  with  a larger 
step-s i ze. 

The  entry  flow  in  a curved  channel  of  constant  centerline  curvature  cr^,  = 
1/ti  has  been  calculated  when  the  channel  height  is  given  by  the  above  expression 
with  S’ | = bn.  For  lip  = 1 or  1.5  accurate  numerical  solutions  to  the  flow  arc 
obtained.  However  with  Ilf  = 2,  the  flow  does  not  separate  but  the  solutions 
oscillate  when  the  step-sizes  are  made  too  small.  For  Ixxuidary  layers  that 
interact  with  the  inviscid  flow,  a similar  lower  limit  on  the  marching  step- 
size  is  encountered  to  avoid  branching  solutions.  These  oscillations  possibly 
indicate  the  slender  channel  equations  are  not  well  posed  as  an  initial  value 
problem  as  multiple  solutions  may  be  possible.  A better  understanding  of  the 
significance  of  these  oscillations  is  needed . 

Rl  I F.RI  NCI1S 

I.  Blottner,  F.G.  , "Numerical  Solutions  of  Slender  Channel  Laminar  Flows,"  to 
be  published. 


12 


FKLE  INTERACTIONS  IN  UNSTEADY  BOUNDARY  LAYERS 


S . N . Brown 

Department  of  Mathematics 
University  College,  London,  U.K. 


ABSTRACT 


The  extension  of  Lighthill's1  theory  of  upstream  influence  to  unsteady' 
supersonic  boundary  layers  was  made  by  Schneider2  who  obtained  upstream- 
propagating  waves  of  small  amplitude  and  high  frequency.  The  Stewartsor.5 - 
Messiter4  triple-deck  theory  can  be  extended  in  a similar  way  and  applied  to 
unsteady  aerodynamics  for  appropriate  orders  of  magnitude  of  these  parameters. 
The  flow  about  the  trailing  edge  of  a supersonic  oscillating  aerofoil  has  been 
discussed5,  and  for  the  subsonic  case  a consistent  structure  obtained5  which 
verifies  the  usual  assumption  of  zero  loading  at  the  trailing  edge  to  remove 
the  non-uniqueness  in  the  potential  flow.  A related  situation  with 
application  to  acoustics  is  the  problem  of  Orszag  and  Crow7  where  there  are 
three  possible  Kutta  conditions.  Daniels8  concludes  that  the  choice  is  not 
determined  tut  presents  a theory  for  high  Strouhal  numbers  which  shows  good 
agreement  with  the  experiments  of  Bechert  and  Pf izenmaier8 . 


REFERENCES 

1.  Lighthill,  M.J.,  "On  boundary  layers  and  upstream  influence.  II.  Supersonic 
flows  without  separation,"  Proc.  Roy.  Soc.  A2 1 7 . pp.  U68-507,  1953. 

2.  Schneider,  W.,  "Upstream  propagation  of  unsteady  disturbances  in 
supersonic  boundary  layers,"  J.  Fluid  Mech.  63.  pp.  ^65-^85,  1 97^ . 

3.  Stewartson,  K. , "On  the  flow  near  the  trailing  edge  of  a flat  plate  II," 
Mathematika  J[6,  pp.  106-121,  1969. 

U.  Messiter,  A.F.,  "Boundary-layer  flow  near  the  trailing  edge  of  a flat 
plate,"  S.I.A.M. , J.  Appl.  Math.  18,  pp.2h1-257,  1970. 

5.  Daniels,  P.G.,  "The  flow  about  the  trailing  edge  of  a supersonic 
oscillating  aerofoil,"  J.  Fluid  Mech.  JL,  pp.  5**1-557,  1975. 

6.  Brown,  S.N.  & Daniels,  P.G.,  "On  the  viscous  flow  about  the  trailing  edge 
of  a rapidly  oscillating  plate,"  J.  Fluid  Mech.  6j,  pp.  7^*3-761 , 1975- 

7.  Orszag,  S.A.  & Crow,  S.C.,  "instability  of  a vortex  sheet  leaving  a serci- 
infinite  plate,"  Studies  in  Appl.  Math.  pp.  167-181,  1970. 

Daniels,  P.G. , "Trai ling-edge  flows,"  I'h.D.  Thesis,  University  of  London, 
197J+  - 


Bechert,  D.  & l’fizenmaier , E.,  "Optimal  compensation  measurements  on  the 
unsteady  exit  condition  at  a nozzle  discharge  edge,"  J.  Fluid  Mech.  7J_, 
pp.123-1^,  1975. 
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THE  TEHEE-DIMENSIONAL  TRIPLE  DECK 


0.  K.  Burggraf 

Tne  Ohio  State  University 
Columbus,  Ohio  43210 


ABSTRACT 


When  combined  with  a viscous  interaction  condition  coupling 
the  pressure  to  the  displacement  thickness , the  boundary-layer 
equations  appear  to  be  an  accurate  model  for  certain  types  of 
separated  flows  at  moderate  Reynolds  numbers.  Examples  of  tne 
success^of  these  methods  are  given  by  Lees  and  Reeves1  and 
Nielsen:  E'or  large  Reynolds  number,  the  boundary  layer  with 

viscous  interaction  subdivides  itself  into  the  substructure  that 
Stewartson'5  has  named  the  triple  uecK.  Numerous  cases  of  two 
dimensional  triple-deck  solutions  have  been^worked  out,  includ- 
ing those  of  Jenson,  burggraf,  and  Rizzetta,  and  of  Rizzetta'3 
for  the  case  of  a compression  ramp  with  a supersonic  mainstream. 

The  theory  is  here  extended  to  three-dimensional  obstacles  with 
supersonic  external  flow.  Applications  to  swept  cylindrical 
obstacles  and  to  ramps  with  triangular  cross-section  are  discussed. 

REFERENCES 

1.  Lees,  L. , and  Reeves,  B. , AIAA  J.  2,  1907-1920  (1964). 
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4.  Jenson,  R. , Burggraf,  0.,  and  Rizzetta,  D.  , Proc.  4th  Int. 

Conf.  on  Numerical  Methods  in  Fluid  Dyn.,  Lecture  Notes  in 
Physics , 35,  (197b). 

p.  Rizzetta,  D.,  Ph.D.  Dissertation,  The  Ohio  State  University 
(1976). 
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AN  INVERSE  TECHNIQUE  FOR  VISCOUS-INVISCID  INTERACTION 


James  E.  Carter 
NASA  Langley  Research  Center 
Hampton,  VA  23665 


ABSTRACT 

The  problem  of  2-D  laminar  and  turbulent  boundary- layer  separation  for 
flows  up  through  the  transonic  speed  range  is  being  studied.  The  boundary- 
layer  equations  for  the  viscous  flow  and  the  small  disturbance  potential 
equation  for  the  inviscid  flow  are  solved  by  novel  techniques  which  are 
specifically  designed  to  compute  separated  flows  with  strong  viscous-inviscid 
interaction.  This  procedure  is  approximate  but  should  provide  reasonably 
accurate  results  for  small  separated  regions,  similar  to  that  found  in  super- 
sonic viscous-inviscid  interactions J Despite  the  notorious  difficulty  in 
numerically  matching  boundary-layer  and  inviscid  solutions,  it  is  thought 
that  the  present  approach,  if  successful,  is  preferable  to  a solution  of  the 
full  Navi er-Stokes  equations  due  to  lower  computer  costs.  However,  it  must 
be  remembered  that  ultimately  the  success  of  either  of  these  approaches 
depends  on  our  capability  to  model  the  turbulent  flow  in  regions  of  large 
adverse  pressure  gradient  including  separation.  This  topic  will  not  be 
discussed  further  here  since  the  present  effort,  so  far,  has  been  directed 
to  the  mathematical/numerical  difficulties  of  iteratively  solving  the  boundary- 
layer  and  inviscid  flow  equations. 

A regular  boundary-layer  solution  is  obtained  at  the  separation  point  by 
prescribing  the  displacement  thickness  as  an  input  condition  to  the  boundary- 
layer  equations;  the  pressure  is  deduced  from  the  resulting  solution.  Catherall 
and  Mangier?  were  the  first  to  use  this  approach;  later,  Carter3  made  a detailed 
numerical  study  of  this  inverse  approach  using  a somewhat  different  formulation. 
In  this  formulation  the  vorticity,  u-component  of  velocity,  and  a perturbation 
stream  function  are  the  dependent  variables.  The  perturbation  stream  function 
is  defined  by  if  = if  - u (y  - 6*)  and  thus  by  requiring  if  -*•  0 as  y -*  <=° 
the  prescribed  displacement  thickness  controls  the  mass  flow  in  the  boundary 
layer.  Recently  a new  formulation  of  the  boundary-layer  equations  for  a 
prescribed  displacement  thickness  has  been  developed  which  is  simpler  since 
it  requires  only  two  dependent  variables,  u and  if.  In  making  comparisons 
with  turbulent  experimental  data,  it  is  simpler  to  work  with  the  u-component 
of  velocity  instead  of  the  vorticity,  since  the  latter  generally  is  not  measured 
and  must  be  deduced  by  differentiating  the  velocity  profile.  The  streamwise 
momentum  equation  is  now  used  in  place  of  the  vorticity  transport  equation 
which  was  used  in  the  earlier  formulation.  The  Crank-Nicholson  finite- 
difference  scheme,  modified  by  the  Reyhner-Fliigge-Lotz  approximation  in 
regions  of  reversed  flow,  is  applied  to  these  equations  with  Newton  linear- 
ization used  to  deduce  the  coupled  linear  algebraic  equations.  The  unknown 
pressure  gradient  is  deduced  implicitly  by  a simple  modification  to  the  Thomas 
algorithm  for  block  tridiagonal  systems.  Alternately,  this  formulation  can 
be  solved  in  the  usual  direct  approach  with  the  edge  velocity,  ue,  prescribed 
instead  of  the  condition  if  = 0 at  the  boundary-layer  edge. 

An  outline  of  this  new  inverse  boundary-1 ayer  procedure  will  be  presented 
along  with  comparisons  to  the  turbulent  boundary-layer  measurements  made  by 
Chu  and  Young  near  the  separation  point.  This  new  procedure  converges  quad- 
ratically  in  3 or  4 iterations  at  each  streamwise  location  for  laminar  flow 
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and  in  less  than  10  iterations  for  turbulent  flow  when  the  maximum  change  in 
either  of  the  dependent  variables  is  required  to  be  less  than  10_5  between  2 
successive  iterations.  A partial  Newton  linearization  on  the  inner  eddy 
viscosity  law  was  found  to  increase  the  convergence  rate  in  the  turbulent 
calculations. 

This  inverse  boundary-layer  technique,  suitably  modified  for  compressibility 
by  the  Stewartson  transformation,  is  currently  being  combined  with  an  inverse 
inviscid  program  developed  by  Carlson. 4 This  inviscid  computation  solves  the 
small  disturbance  potential  equation  for  transonic  flow  for  a prescribed 
surface  pressure  distribution  which  is  provided  by  the  inverse  boundary-layer 
solution.  The  inviscid  solution  then  gives  the  slope  of  displacement  body 
which  is  integrated  and  combined  with  the  actual  surface  ordinates  to  yield 
the  displacement  thickness  for  the  input  condition  to  next  boundary-layer 
computation.  The  form  of  this  inverse-inverse  iterative  procedure  is 
identical  to  that  presented  by  Carter  and  Wornom^  which  was  used  successfully 
to  solve  a low-speed  laminar  separation  problem.  In  that  computation  it  was 
found  that  no  numerical  smoothing  was  required  as  is  typically  used  in  matching 
boundary-layer  and  inviscid  flows.  This  result  occurs  because  classical 
matching  requires  numerical  differentiation,  whereas  the  present  procedure 
uses  numerical  integration  of  the  computed  inviscid  and  boundary-layer  results. 

At  the  present  time,  Carlson's  program,  now  modified  to  include  solid  wind- 
tunnel  walls,  is  being  used  to  compute  the  inviscid,  transonic  flow  over  the 
bump  on  the  lower  tunnel  wall  which  was  measured  by  Alber,  et  al.6  The  experi- 
mental pressure  is  used  as  the  input  condition  and  it  is  found  that  the  computed 
displacement  thickness  agrees  fairly  well  with  that  measured  by  Alber,  et  al. 
Further  calculations  are  required,  though,  as  it  has  been  observed  that  the 
deduced  displacement  thickness  on  the  bump  on  the  lower  tunnel  wall  is  very 
sensitive  to  the  slope  of  the  upper  tunnel  wall.  Hence,  at  the  present  time 
it  appears  that  it  will  be  necessary  to  include,  in  the  usual  fashion,  the 
boundary-layer  displacement  thickness  on  the  upper  tunnel  wall.  For  the  work- 
shop the  current  status  of  the  application  of  this  inverse  interaction  technique 
will  be  presented. 
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SEPARATED  FLOWS  AND  THEIR  REPRESENTATION  BY  BOUNDARY-LAYER  EQUATIONS* 

Tuncer  Cebeci 

California  State  University  at  Long  Beach 
Long  Beach,  California,  90815,  U.S.A. 

ABSTRACT 


The  paper  describes  two  separate  approaches  for  calculating  laminar  and 
turbulent  boundary  layers  with  positive  and/or  negative  wall  shear.  The  first 
approach,  which  we  call  the  "nonlinear  eigenvalue  method,"  treats  the  pressure 
as  an  eigenvalue  parameter  and  works  with  either  physical  or  transformed  variables. 
The  second  approach  treats  the  pressure  as  an  unknown  function  and  seeks  the  solu- 
tion of  the  governing  equations  in  physical  variables.  This  approach,  which  we 
call  the  "Mechul"  function  scheme,  was  used  earlier  for  the  Falkner-Skan  equation 
by  Cebeci  and  Keller.  Our  studies  show  that,  while  the  first  approach  only  works 
for  flows  with  positive  wall  shear,  the  second  approach  works  extremely  well  for 
flows  with  either  positive  or  negative  wall  shear.  This  abstract  concentrates  on 
the  second  approach. 

The  Mechul  method  is  relevant  to  the  prediction  of  separation  and  the  calcu- 
lation of  the  flow  properties  within  separated  flow  regions.  It  has  been  developed 
specifically  to  allow  the  determination  of  flow  properties  in  the  vicinity  of  the 
suction  side  of  airfoils.  Boundary- 1 ayer  equations  are  solved  and  are  appropriate 
upstream  of  separation  and  more  approximate  in  the  separated  flow.  As  is  demon- 
strated, however,  when  used  in  the  manner  described  in  this  paper  and  outlined 
below,  their  solution  can  lead  to  results  which  are  sufficiently  precise  for  design 
purposes  and  require  very  short  computer  times  and  little  storage.  The  method  also 
allows  the  testing  of  various  turbulence  models  for  separated  flows. 

A brief  description  of  the  Mechul  function  method  for  specified  displacement 
thickness  follows.  With  the  concept  of  eddy  viscosity,  em,  and  stream  function, 

0,  the  boundary-layer  equations  for  two-dimensional  laminar  and  turbulent  flows 


can  be  written  in  the  form, 

(b*")  ~ Px  = If  H) 

where  primes  denote  differentiation  with  respect  to  y and  b = 1 + em/v.  Eq.  (1) 
is  subject  to  the  boundary  conditions 

y = 0 ; ip  = i (/ 1 = 0 (2a) 

y = ye;  = ue  , = ug[ye  - 6*(x)]  (2b) 


The  pressure  p is  treated  as  an  unknown  and  equation  (1)  is  written  as  the 
first-order  system, 


<(/'  = u 


(3a) 


u'  = v (3b) 


*This  work  was  supported  by  the  Office  of  Naval 
N0001 4-74-A-0203-0001  . 
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P'  = o (3c) 

(b*)'  -px  -uf-.ft  (3d) 

Similarly,  eqs.  (2)  become 

y = 0;  ik  = u = 0 (4a) 

y = ye;  u = ug  , <Le  = ue[ye  - 5*(x)]  (4b) 


The  solution  of  the  system  (3)  and  (4)  is  obtained  with  Keller's  two-point 
finite-difference  method.  In  the  region  of  negative  longitudinal  velocity,  the 
FLARE  approximation  discussed  by  P.  G.  Williams  was  used. 

Several  flows,  which  allow  the  displacement  thickness  distribution  to  be 
specified,  have  been  calculated  with  this  method.  These  include  the  laminar 
separated  flows  of  P.  G.  Williams,  J.  E.  Carter,  and  R.  Briley;  In  the  last  case, 
it  was  possible  to  make  comparison  with  solutions  of  the  Navier-Stokes  equations 
and  the  agreement  is  shown  to  be  excellent.  The  experimentally  investigated  flow 
of  R.  Simpson,  et  al .was  also  calculated  and  a comparison  of  the  1 aser-annemometer 
measurements  with  the  present  results  is  presented  and  shown  to  be  within  the 
limits  of  experimental  uncertainty. 
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STUDY  ON  UNSTEADY  TRA I L I NG-EDGE  STALL 


H.K.  Cheng 

University  of  Southern  California 
Los  Angeles,  California  90007  USA 


ABSTRACT 


The  asymptotic  theory  of  i nv i sc i d-vi scous  interaction  near  a sharp  trailing 
edge  is  studied  from  the  viewpoint  of  aquatic  animal  propulsion.  Development  of 
thrust  on  the  caudal  fin  (as  well  as  on  a bird's  wing)  depends  largely  on  the 
(unsteady)  Kutta  condition,  aside  from  an  unseparated  smooth  flow  around  the 
leading  edge.  The  departure  from  the  ideal  Kutta  condition  (of  vanishing  pressure 
jump)  at  the  edge  at  high  and  low  (reduced)  osc i 1 1 at i ng  f requenc i es  is  of  great  interest. 

The  theory  for  the  case  corresponding  to  a flat  plate  at  incidence  od  * has 
been  set  forth  by  Brown  and  Stewartsonl,  adopting  the  basic  three-deck  structure  of  fef  1 
around  the  trailing  edge  (T.  E.)  The  theory,  as  a self-consistent 

asymptotic  analysis,  hold> as  long  as 


= C’(&  , 


where  (f  = ' . An  extension  to  the  case  of  a rapidly  oscillating  plate 

has  been  made  by  Brown  and  Daniels^.  To  match  the  three-deck  structure^a  Stokes 
wall  layer  and  two  fore  decks  are  introduced  in  Ref.  3-  Whereas  interesting 
viscous  correction  to  the  circulation  has  been  determined  (estimated)  as  in  Ref.  1 
for  the  steady  case,  requirements  on  * and  the  frequency  for  the  model  appear 
to  be  rather  severe.  Namely  , 


ct*  - G i€  J 

SS^  = 0(«!)  J 


We  note  that, in  most  problems  of  aero  / hydromechanics,  cd  must  be  greater  than 
the  main-deck  thickness  W l.  ~ Ofe<*^and  that  the  range  £ =Q(l)  is  of  great 
practical  interest.  The  work  in  Ref.  3 nevertheless  brings  out  much  insight  to 
the  problems  and  furnishes  a valuable  guide  to  the  analysis  for  other  domains  of 
^ * and  5 . 

The  study  on  the  unsteady  trailing  flow  to  be  carried  out  will  focus  on 
ranges  larger  amplitude  and  lower  frequency,  especially  domains  with 

| ) 

O - 0(1)  . j (3) 

One  of  these  domains,  namely  , 

a"  = O te*)  ) 

appears  tractable,  for  which  a structure  similar  to,  but  slightly  simplier  than, 
Brown  and  Dan i e 1 s ' ^appears  to  be  applicable. 


It  is  interesting  to  examine  the  high  and  low  frequency  effects  as  a conse- 
quence of  shorf^long  wave  lengths  (compared  to  w:ng  chord).  Fcjr^high  frequency 
( 5'  >>  / ) , the  influence  of  the  periodic  tva|<£  vorticity  averages^out  (to  a weaker 
level),  and  the  velocity  and  pressure  on  the  wing  (lit of  the  T.  E.)  is  dominated  by 
the  apparent-mass  and  the  quasi-steady  effects.  This  is  found  to  be  in- 

deed tne  case  in  Ref.  3,in  which  the  wake-vort i c i ty  effect  does  not  enter  into 
the  analysis,at  least  for  the  fore  decks.  We  expect  a more  pronounced  periodic 
wake  influence  for  the  case  under  condition  (4).  The  presentation  will  describe 
a preliminary  examination  of  the  asymptotic  structure  and  scaling  problems. 
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BOUNDARY-LAYER  CALCULATIONS 


FOR 

A CIRCULAR  CYLINDER 
S.C.R.  Dennis 

The  University  of  Western  Ontario 
London,  Ontario,  Canada 

ABSTRACT 

In  order  to  solve  the  two-dimensional  laminar  boundary- layer  equations  it 
is  necessary  to  specify  the  external  pressure  distribution  at  the  outer  edge  of 
the  boundary  layer  or  some  equivalent  information.  It  is  clear  from  the 
tendency  of  the  surface  pressure  distributions  obtained  from  numerical  solutions 
of  the  Navier-Stokes  equations  as  the  Reynolds  number  increases  that  the 
external  pressure  distribution  given  by  potential  flow  is  not  an  appropriate 
assumption.  In  the  present  note  some  results  are  presented  which  have  been 
calculated  by  integrating  the  boundary-layer  equations  for  a circular  cylinder 
using  pressure  distributions  obtained  from  previous  numerical  solutions  of  the 
Navier-Stokes  equations.  The  object  is  to  check  the  results  for  consistency 
with  the  trend  of  the  numerical  solutions  for  increasing  Reynolds  number.  The 
work  is  a first  step  in  an  attempt  to  devise  an  interaction  scheme  between  the 
boundary-layer  equations  and  a model  of  the  external  flow  obtained  from  the 
Navier-Stokes  equations. 

The  integrations  of  the  boundary-layer  equations  have  been  carried  out  in 
three  separate  cases  utilizing  pressure  distributions  taken  from  the  numerical 
solutions  of  Dennis  & Chang^  at  R = 40,  70  and  100  respectively.  Here 
R = 2U00a/v,  where  a is  the  radius  of  the  cylinder,  Uco  the  stream  velocity  arid  v 
tire  coefficient  of  kinematic  viscosity.  A calculation  of  a similar  nature  has 
previously  been  given  by  Dimopoulos  & Hanratty^  for  the  single  case  R = 40 
utilizing  the  pressure  distribution  calculated  from  a numerical  solution  of  the 
Navier-Stokes  equations  by  Kawaguti  & Jain3.  The  present  work  is,  however,  able 
to  estaoiish  some  notion  of  a trend  of  the  results  with  the  changing  external 
pressure  distribution  as  R increases.  There  is  found  to  be  a good  measure  of 
consistency  between  the  boundary-layer  calculations  and  the  numerical  solutions 
of  the  Navier-Stokes  equations  of  Dennis  & Chang. 

Some  results  are  given  for  the  variation  of  local  skin  friction  over  the 
surface  of  the  cylinder  and  the  variation  of  velocity  throughout  the  boundary 
layer.  There  is  a striking  illustration  of  the  breakdown  of  the  boundary-layer 
equations  near  the  point  of  separation.  The  consistency  of  the  boundary-layer 
calculations  with  the  numerical  solutions  of  the  full  Navier-Stokes  equations 
near  the  front  stagnation  point  is  illustrated  by  the  following  results.  The 
pressure  coefficient  at  the  front  stagnation  point  is  given  by  boundary- laver 
theory  in  the  form 

(P0  ' Pj/K  ^ 1 + a/R, 

where  is  the  density  and  p , pn  the  pressures  at  the  I ronL  stagnation  point 
and  at  large  distances,  respectively.  The  estimate  a = 6.35  is  obtained  from  all 
three  boundary-layer  solutions,  which  is  not  inconsistent  with  the  trend  of  the 
estimates  i = j.18,  5.76,  5.95,  6.00  obtained^respectively  at  R = 20,  40,  70  and 
100  by  Dennis  \ Chang.  Finally,  if  cf  = xj  is  the  local  coefficient  of 
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skin  friction  at  the  surface  of  the  cylinder,  where  x is  the  local  shear 
stress,  the  slope  of  the  curve  of  cf  at  the  front  Stagnation  point  with 
respect  to  angular  distance  round  the  cylinder  comes  out  to  be  6.97  from  all 
three  boundary-layer  solutions.  This  is  quite  consistent  with  the  estimates 
6.59,  6.70  and  6.89  of  this  quantity  given  at  R = 40,  70  and  100  by  Dennis  & 
Chang. 
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FLOW  IN  A CURVED  PIPE  AT  HIGH  DEAN  NUMBER 


S.  C.  R.  Dennis 

University  of  Western  Ontario,  London,  Canada 

and 

N.  Riley 

University  of  East  Anglia,  Norwich,  U.K. 
ABSTRACT 


The  flow  down  a circular  pipe  of  radius  a,  coiled  in  a circle  of  radius  L, 
under  the  action  of  a constant  axial  pressure  gradient  C,  is  considered.  The 
first  successful  analysis  of  such  a flow  was  presented  by  Deanl>2  for  the  case 
of  small  values  of  the  parameter  D = (Ca2/p) (2a3/v2L) : , where  p,  v represent 
the  viscosity  and  kinematic  viscosity  of  the  fluid  respectively,  and  vanishingly 
small  a/L.  Our  concern  is  with  the  corresponding  flow  for  D >>  1,  If, 
following  Collins  and  Dennis3,  we  write  the  stream  function  for  the  (secondary) 
velocity  components  in  a cross-flow  plane  as  T = vij)  and  the  axial  velocity  as 
W = (v2L/2a3)5  w then  to  leading  order  an  inviscid  'core'  solution  is  obtained 
for  which  ip  = D/1r  sin  0/w^,(r  cos  0) , w = D^3 wc(r  cos  0).  The  core  velocity 
w (r  cos  0)  is  not  determined  at  this  stage.  To  determine  wc  and  satisfy 
t§e  conditions  at  r = 1 it  is  necessary  to  introduce  a boundary  layer  close  to 
r = 1.  Classical  arguments  indicate  that  this  is  of  thickness  0(aD-l/3),  which 
is  in  accord  with  numerical  results3  for  large  but  finite  D.  The  core  velocity 
wc  which  provides  an  outer  boundary  condition  for  the  axial  velocity  within  the 
boundary  layer  is  itself  determined  by  ensuring  continuity  of  mass  flow  between 
the  core  and  boundary  layer.  In  this  sense  the  problem  is  an  interaction 
problem.  By  exploiting  the  mass  flow  continuity  condition,  wc  may  be  formally 
eliminated  from  the  problem,  except  that  it  must  be  given  a prescribed  value, 
wc  say,  at  the  outside  bend  0=0  (see  figure).  The  asymptotic  structure 
close  to  0 = it  is  known,  and  in  principle  the  solution  can  be  advanced  step  by 
step  from  0 = 0;  wc  is  determined  at  each  step.  The  correct  value  of  wc  is 
that  for  which  the  numerical  solution  approaches  the  asymptotic  form  smooFhly 
as  0 -*■  v . Difficulties  have  been  experienced  with  the  above  scheme,  since 
although  separation  of  the  secondary  cross-flow  has  not  been  encountered  a 
region  of  reversed  flow  close  to  the  outer  edge  of  the  boundary  layer  has. 

This  leads  to  a breakdown  of  the  forward  integration  scheme,  and  an  attempt  has 
been  made  to  treat  the  problem,  numerically,  as  a boundary-value  problem.  The 
results  which  have  been  obtained  are  encouraging  but,  to  date,  incomplete. 
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A CONTROLLING  FACTOR  IN  THE  ESTABLISHMENT  OF  THE 


SIZE  OF  THE  SEPARATION  REGION 

W.  L.  Hankey  and  J.  S.  Shang 

Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  AF  Base,  Ohio  45433 

ABSTRACT 


Solutions  of  the  complete  Navier  Stokes  equations  have  been  obtained  for 
shock-laminar  boundary  layer  interactions  by  using  MacCormacks  finite  differ- 
ence scheme.  Turbulent  interactions  also  have  been  computed  with  this  program 
by  substitution  cf  an  eddy  viscosity  into  the  equations  and  tireatj.ug  piuperLies 
as  the  time  averaged  mean  values.-*  It  was  found  that  a relaxation  of  the  eddy 
viscosity  with  a "time  constant"  of  A=106  was  necessary  to  predict  the  experi- 
mentally observed  separation  length.  To  further  examine  the  details  of  this 
interaction,  the  total  pressure  profiles  were  compared  before  and  after  the 
separation  region  (Fig.  1).  In  the  outer  inviscid  region  large  shock  losses 
are  observed,  however,  an  increase  in  total  pressure  along  the  wall  is  noted. 
Since  static  pressure  increases  through  the  interaction  the  total  must  also 
rise  in  that  vicinity.  To  answer  the  question,  "under  what  conditions  can  total 
pressure  increase  in  a viscous  region?",  the  energy  equation  was  examined.  In- 
tegration of  the  energy  equation  along  an  adiabatic,  steady,  boundary  layer 
streamline  produces  the  following  relationship: 

P x9 

02  r 2 

In  = \ T„  dx 

This  integral  equation  shows  that  total  pressure  can  be  changed  according  to  the 
sign  of  Ty  A local  decrease  in  entropy  is  possible,  however  the  net  total 
change  in  entropy  must  of  course  increase  (in  accord  with  the  second  law  of 
thermodynamics) . 

Much  can  be  deduced  from  the  integral  relationship  concerning  the  shape  of 
the  velocity  profile  in  the  vicinity  of  the  recirculation  region.  For  example, 
the  flow  just  above  the  dividing  streamline  experiences  an  increase  in  total 
pressure  through  the  interaction.  Therefore  Ty  must  be  positive  in  some  region 
above  the  dividing  streamline.  Also  it  is  obvious  that  some  length  (x2  - x^) 
is  necessary  to  provide  sufficient  time  for  the  viscous  forces  to  act  to  in- 
crease Po . 

Examination  of  similar  flows  confirms  this  supposition.  Figure  2 snows  the 
locus  of  the  inflection  points  for  similar  reversed  flows  which  divide  the  re- 
gions of  negative  and  positive  Ty  Note  that  a small  region  does  exist  above 
the  dividing  streamline  for  which  Ty>0.  Flow  in  this  region  will  experience  an 
increase  in  total  pressure,  while  flow  above  will  only  undergo  a decrease  in 
total  pressure. 


By  assuming  locally  similar  flows  it  is  possible  to  estimate  the  extent  of 
laminar  separation  by  employing  the  previously  derived  integral  relationship  for 
total  pressure  change.  By  using  a mean  value  for  the  integrand  an  approximation 
for  the  length  of  separation  is  derived  which  is  similar  to  the  correlation  de- 
duced by  Ball^  from  six  different  experimental  investigations. 

Also  the  results  of  the  turbulent  relaxation  model  may  be  analyzed  using 
this  integral  equation.  In  reference  1 the  separation  length  was  observed  to 
increase  dramatically  from  the  equilibrium  value  when  the  eddy  viscosity  was 
frozen  at  the  upstream  value  thereby  decreasing  Tv  near  the  dividing  streamline. 
One  concludes  that  the  length  of  separation  is  extremely  sensitive  to  Ty  only  in 
the  vicinity  of  the  dividing  streamline  and  therefore  the  value  of  the  eddy 
viscosity  in  other  regions  is  not  nearly  as  important.  For  example,  Shang  found 
the  separation  length  to  vary  only  slightly  when  the  eddy  viscosity  in  the  sepa- 
ration bubble  was  changed  between  zero  and  the  frozen  value.  Therefore  it  can 
be  concluded  that  what  goes  on  in  the  recirculation  zone  is  unimportant  in  es- 
tablishing the  length  of  the  separation  zone  and  analyses  of  only  tecirculacion 
then  cannot  lead  to  any  strong  condition  (this  has  been  perplexing  for  years). 
However  an  analysis  in  the  region  immediately  above  the  recirculation  region 
does  produce  a useful  result. 

In  summary,  a simple  integral  relation  produces  a condition  which  implies 
that  Ty  must  be  positive  above  the  dividing  streamline  in  order  to  permit  a de- 
crease in  entropy  (or  an  increase  in  total  pressure)  in  the  lower  portion  of 
the  boundary  layer.  This  relationship  is  felt  to  be  the  single'most  controlling 
factor  in  establishing  the  size  of  the  separation  region.  Only  analyses  which 
properly  account  for  this  phenomenon  can  produce  accurate  results. 
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RECENT  EXPERIMENTAL  STUDIES  OF  THE  MECHANISM  OF  FLOW  SEPARATION  IN 
HIGH  MACH  NUMBER,  HIGH  REYNOLDS  NUMBER  FLOWS 


Michael  S.  Holden 

Calspan  Corporation 
Buffalo,  New  York  1*4221 


ABSTRACT 

A discussion  is  given  of  the  results  of  recent  experimental 
studies  which  have  important  implications  on  the  modeling  of  lami- 
nar and  turbulent  separated  flows,  and  present  new  problems  in  this 
area  of  research.  The  results  of  detailed  experimental  studies  of 
both  wedge-  and  shock-induced  regions  of  viscous  interaction  and 
boundary  layer  separation  are  discussed  in  terms  of  their  implica- 
tions to  the  theoretical  modeling  of  these  flows.  The  studies  of 
laminar  separation  induced  "internally"  at  the  junction  between  a 
flat  plate  and  a wedge,  and  "externally"  with  an  oblique  shock, 
demonstrated  that  at  Mach  numbers  from  12  to  18  the  laminar  bound- 
ary laye . can  be  more  easily  separated  by  an  externally  generated 
snock.  This  result  is  in  contrast  to  earlier  findings  at  lower 
Mach  numoers,  and  suggests  that  modeling  these  flows  with  a simple 
"free  interaction"  theory  may  be  inadequate.  Studies  of  shock- 
and  wedge-induced  interaction  regions  at  Mach  numbers  from  6 to  13 
and  Reynolds  numbers  up  to  200  x 106  suggest  that  the  pressure  rise 
to  separation  is  independent  of  the  mechanism  inducing  separation. 
These  studies  suggest  that  boundary  layer  separation  is  controlled 
primarily  oy  an  interaction  wnich  occurs  at  the  base  of  the  bound- 
ary layer,  thus  describing  the  separation  mechanism  in  terms  of  the 
interaction  between  the  growtn  of  a turbulent  boundary  layer  and 
an  external  inviscid  flow  is  inaccurate  at  high  Mach  numbers.  The 
length  of  these  turbulent  interaction  regions  are  comparable  with 
ooundary  layer  thickness,  thus  it  is  questionable  whether  first 
order  boundary  layer  theory  is  valid  in  these  flows.  Studies  of 
steady  and  unsteady  separated  regions  over  concave  nose  shapes 
generated  by  non-similar  ablation  are  discussed  briefly.  Finally, 
we  uiscuss  some  novel  studies  wiiere  extensive  regions  of  unsteady 
separated  flows  were  developed  when  a minute  particle  launched 
tnrough  tne  stagnation  region  of  a blunt  oody  pierced  the  Low 
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UPSTREAM  INFLUENCE  IN  INTERACTING 
NONSEPARATING  TURBULENT  BOUNDARY  LAYERS 
G.  R.  Inger 

Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Va.,  24061,  U.S.A. 

ABSTRACT 


The  important  feature  of  upstream  influence  in  all  viscous-inviscid 
interaction  problems  is  known  to  be  primarily  governed  by  the  viscous 
behavior  of  the  interactive  disturbances  within  a thin  "frictional 
sub-layer"  near  the  body  surface.!  Existing  analytical  treatments  of 
this  region  assume  that  it  lies  within  the  linear  portion  of  the  incoming 
flow  velocity  profile.  Although  valid  for  laminar  flows,  this  assumption 
fails  badly  for  turbulent  flows  at  the  higher  Reynolds  numbers  Re^  = 106-1C)!0 
typifying  full  scale  flight  vehicles.  In  the  present  paper  we  present  the 
formal  extension  and  complete  solution  for  the  turbulent  non-separating 
case  wherein  the  velocity  profile  is  not  necessarily  linear  (i.  e.,  the 
frictional  sublayer  may  penetrate  to  an  arbitrary  degree  into  the  loga- 
rithmic part  of  the  profile.  ) 
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Freely  Interacting  Transonic  Boundary  Layers 


★ 


A.  Kluwick 

Technische  Universitat  Wien 
and 

R.  j.  Bodonyi 

Indiana  University  - Purdue  University  at  Indianapolis 


The  inner  layer  equations  of  the  triple  deck  structure  for  transonic 
free  interactions! , both  compressive  and  expansive,  have  been  studied 
numerically  and  analytically.  The  only  restriction  placed  on  the  external 
flow  is  that  it  remains  supersonic  throughout. 

In  the  case  of  compressive  interactions  characterized  by  values  of 
the  transonic  interaction  parameter  Ko  < 2.3  the  assumption  of  supersonic 
external  flows  breaks  down  before  separation  occurs.  For  Ko  > 2.3,  how- 
ever, the  supersonic  flow  region  extends  into  the  region  downstream  oi  ^ 
the  separation  point.  Calculations  based  on  the  Reyhner  and  Flugge-Lotz^ 
approximation  in  the  reversed  flow  reqion  indicate  that  for  sufficiently 
larqe  values  of  Ko  the  pressure  distributions  approach  regions  of  constant 
pressure  which  is  similar  to  the  supersonic  flow  case^. 

The  results  also  show  that  the  ultimate  upstream  decay  is  identical 
to  that  of  the  supersonic  free  interaction  if  Kq  0.  For  expansive  free 
interactions  with  Ko  = 9,  however,  the  upstream  decay  is  algebraic  rather 
than  exponential . 

Finally,  the  asymptotic  structure  of  the  flow  field  near  the 
pressure  singularity  in  transonic  expansive  flows  has  been  determined,  and 
the  structure  agrees  well  with  the  numerical  computations  near  the  sing- 
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CALCULATION  OF  AXISYMMETRIC  SEPARATED 
TURBULENT  BOUNDARY  LAYERS 

Gary  D.  Kuhn 

Nielsen  Engineering  & Research,  Inc. 

Mountain  View,  California  94043  U.S.A. 

A method  is  being  developed  to  analyze  compressible  viscous 
attached  and  separated  flow  fields  for  axisymmetric  body 
geom  :ries  used  for  isolated  nozzle  afterbody  models.  The  work 
includes  coupling  of  the  integral  boundary-layer  method  with  a 
finite-difference,  inviscid-f low  program  using  the  boundary-layer 
displacement  thickness  to  modify  the  shape  of  the  body. 

The  derivation  of  the  governing  equations  for  the 
boundary-layer  flow  has  been  described  in  reference  1.  That 
derivation  has  been  modified' to  use  the  integral  method  described 
in  reference  2.  The  inviscid  flow  is  calculated  using  the  South 
and  Jameson  program  described  in  reference  3.  That  program  is 
also  being  modified  to  produce  a finite-difference  solution  of 
the  full  potential  equation  in  conservative  form. 

The  boundary-layer  method  uses  the  law-of-the-wall  and  the 
law-of-the-wake  velocity  profile  formulation  with  a laminar  sub- 
layer included.  This  formulation  is  capable  of  calculating  both 
laminar  and  turbulent  boundary  layers  as  well  as  separated 
boundary  layers.  For  attached  flow,  the  boundary-layer  edge 
velocity,  U , is  prescribed  and  the  boundary-layer  thickness, 

6,  and  the  friction  velocity,  u , are  calculated.  For  separated 
flows,  th.e  displacement  thickness,  6*,  is  prescribed  and  the 
thickness,  6,  and  the  edge  velocity,  UQ , are  calculated. 

The  axisymmetric,  conservative,  finite-difference  inviscid 
flow  program  will  solve  the  equation  for  conservation  of  mass 
employing  artificial  viscosity  to  assure  upwind  differencing  in 
supersonic  zones.  The  particular  form  of  the  artificial  viscosity 
and  the  iterative  relaxation  procedure  to  be  used  are  currently 
under  development. 


Several  kinds  of  iterative  procedures  are  being  developed 
for  calculating  the  interaction  of  the  viscous  and  inviscid  flows. 
In  one  approach,  the  displacement 'thickness  is  prescribed  on  the 
boundary  layer,  and  the  free-stream  velocity  distribution  which 
results  from  the  boundary-layer  calculation  is  compared  with  that 
from  the  inviscid  flow  calculation  on  the  augmented  body.  The 
displacement  thickness  is  then  modified  by  trial  and  error  until 
the  two  pressure  distributions  agree.  Other  approaches  are  being 
studied  which  are  expected  to  lead  to  a calculation  method  which 
will  iterate  without  external  assistance. 

The  work  described  in  this  abstract  is  being  sponsored  by 
the  Arnold  Engineering  Development  Center,  Air  Force  Systems 
Command,  USAF. 
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AN  EFFICIENT  NUMERICAL  METHOD  FOR  SOLVING  THE  TIME-DEPENDENT 
COMPRESSIBLE  NAVIER-STOKES  EQUATIONS  AT  HIGH  REYNOLDS  NUMBER 


Robert  W.  MacCormack 

Ames  Research  Center,  NASA 
Moffett  Field,  Calif.  94035  USA 


ABSTRACT 

The  Navier-Stokes  equations  adequately  describe  aerodynamic  flows  at  stan- 
dard atmospheric  temperatures  and  pressures.  If  we  could  efficiently  solve  these 
equations,  there  would  be  no  need  for  experimental  tests  to  design  flight  vehicles 
or  other  aerodynamic  devices.  Unfortunately,  analytic  or  closed-form  solutions 
to  these  equations  exist  for  only  a few  simple  flow  problems.  During  the  past 
decade,  the  computer  has  been  used  to  generate  many  new  solutions.  However,  with 
existing  numerical  methods1  and  computer  resources,  these  solutions  have  been 
restricted  to  low  Reynolds  number  or  two-dimensional  flows.  This  paper  describes 
a new  numerical  method"  that  has  been  used  to  drastically  reduce  the  computation 
time  required  to  solve  the  Navier-Stokes  equations  at  flight  Reynolds  numbers. 
Though  flows  past  complete  aircraft  configurations  are  still  beyond  our  reach, 
the  new  method  makes  it  possible  and  practical  to  calculate  many  important  three- 
dimensional,  high  Reynolds  number  flow  fields  on  today's  computers. 

The  unsteady  Navier-Stokes  equations  are  difficult  to  solve  because  at  high 
Reynolds  numbers  the  magnitude  of  the  inertial  forces  described  by  the  hyperbolic 
terms  of  the  equations  are  much  larger  than  the  viscous  forces  described  by  the 
parabolic  terms.  Conventional  numerical  methods,  whether  explicit  or  implicit, 
are  severely  restricted  to  small  time  steps  by  the  stability  or  accuracy  require- 
ments imposed  by  the  hyperbolic  terms;  hence,  many  time  steps  may  be  required 
before  the  viscous  effects  can  be  determined.  The  present  approach  time-splits 
the  equations  into  a hyperbolic  part  and  a parabolic  part,  solves  the  hyperbolic 
part  by  using  a new  explicit  numerical  method  based  on  characteristics  theory, 
and  solves  the  parabolic  part  by  using  a new  efficient  implicit  parabolic  method. 
Both  methods  are  fully  conservative  and  stable  with  time-steps  orders  of  magni- 
tude larger  than  those  allowed  by  CFL  (Courant,  Friedrich,  and  Lewy)  stability 
criteria. 
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VISCOUS  INTERACTION  AND  BOUNDARY- LAYER  SEPARATION 


H.  McDonald*  and  W.  R.  Briley** 
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ABSTRACT 

Given  the  requirements  of  having  to  predict  the  details  of  the  viscous 
interactions  and  separated  flows  over  a very  wide  range  of  flow  conditions, 
the  present  authors  have  concentrated  on  making  predictions  by  computing 
efficient  and  direct  numerical  solutions  to  the  governing  partial  differential 
equations  of  motion.  Aspects  of  relatively  low  Reynolds  number  laminar 
incompressible  separation  bubbles  have  been  investigated  using  this  approach 
by  Brileyl.  Of  note  in  this  work  is  the  observation  of  considerable  upstream 
influence  of  the  separation  bubble  at  these  relatively  low  Reynolds  numbers 
and  the  apparent  absence  of  any  separation  singularity.  Later  the  incompressible 
numerical  . ieme  was  further  developed  and  improved  and  higher  Reynolds  number 
: •.minar  transitional  and  turbulent  incompressible  separation  bubbles  computed 
in  '•  time  dependent  manner  by  Briley  and  McDonald^.  Both  the  full  Navier- 
Gtc'fc.:  equations  and  the  boundary  layer  equations  were  used  and  compared.  Of 
iii  te  m this  work  was  the  use  of  the  time  dependent  approach  to  allow  the 
interaction  of  the  shear  flow  region  with  the  outer  inviscid  flow  and  the 
a.  parent  absence  of  the  separation  singularity  when  this  interaction  is  allowed 
for  with  the  boundary  layer  equations.  Further  the  highly  implicit  and  non- 
linear nature  of  the  turbulence  model  turned  out  to  be  amenable  to  treatment 
within  the  computational  framework  and  the  predicted  results  were  in  good  agree- 
ment with  data.  Mure  recent  studies  have  been  concerned  with  the  compressible 
problem  and  here  both  the  steady  boundary  layer  equations  and  the  time  dependent 
Navier-Stokes  equations  arc  being  used.  Levy,  Shamroth  and  Me Donald ^ discuss 
results  obtained  in  both  laminar  and  turbulent  shock  wave  boundary  layer  inter- 
actions containing  small  regions,  of  flow  separation  in  supersonic  flow  using 
iteratively  a spatially  marching  boundary  layer  scheme.  A companion  study  on 
t sc  .vine  sr  LIem  uses  he  fully  implicit  scheme  for  solving  the  multi -dimensional 
Navior-oU  ■■■.:■  ' eqi.a*  : n evei  ined  by  Briley  and  McDonald*4.  Of  particular  note 
in  v A v r r-  tut  very  modest  computer  run  times  attainable  with  the  fully 
‘ Ja . • • -■  ' < .seme,  and  this  should  enable  future  extensive  and 

« . t.  ; .:  rob lem  using  the  full  equations  without  exorbitant 

'•  ‘ " expen.  . lastly  some  preliminary  work  which  should  eventually 

’ • m t utiun  of  tiie  flow  around  an  arbitrary  airfoil  using  the 

r>  i .re  of  Briley  and  Me  Dona  Id1*  will  be  discussed. 
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TURBULENT  INTERACT  IONS  AT  TRAILING  EDGES 


R.  E.  Melnik 
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ABSTRACT 

We  are  considering  the  effect  of  turbulent  boundary  layers  on  the  lifting 
characteristics  of  airfoils  at  high  Reynolds  numbers.  The  most  important  effect 
arises  from  the  local  boundary-layer  behavior  near  the  trailing  edge,  which 
directly  effects  the  Kutta  condition  and  leads  to  significant  reductions  in  lift. 
Existing  methods  for  predicting  the  boundary-layer  effects  on  lift  are  based  on 
the  displacement  surface  corrections  of  a second-order  boundary- layer  theory. 

> The  displacement  surface  is  assumed  to  be  a streamline  and  the  corrections  to  the 
inviscid  flow  are  determined  by  standard  procedures.  Unfortunately,  the  displace- 
ment thickness  distribution  determined  from  the  solution  of  the  boundary-layer 
equations  develops  singularities  at  trailing  edges.  The  singularities  are  caused 
by  the  discontinuous  jump  in  the  surface  boundary  conditions  across  the  trailing 
edge  and  by  a singularity  in  the  pressure  gradient  of  the  inviscid  solution  at  the 
trailing  edge.  Because  of  these  singularities  a viscous  Kutta  condition  cannot  be 
imposed  and  the  lift  coefficient  cannot  be  determined  by  standard  second-order  theory. 

To  remedy  this  situation,  we  have  developed  an  "inner  solution"  for  the 
strong  interaction  region  at  trailing  edges.  In  the  present  discussion,  we  will 
outline  the  main  feature"  of  the  interaction  theory  and  will  describe  recent 
developments  leading  .o  a complete  solution  for  the  lift  decrement  due  to  turbu- 
lent boundary  layers. 

The  interaction  theory  of  Ref.  1 is  based  on  a large  Reynolds  number, 
asymptotic  expansion  of  the  full  (and  unclosed)  Reynolds  equations  of  turbulent 
flow.  The  theory  leads  to  a three-layer  description  of  the  flow  near  the  trail- 
ing edge  consisting  of:  1)  an  outer,  inviscid,  rotational  flow  region  encompass- 
ing most  of  the  boundary  layer,  2)  a thin,  conventional,  constant  stress,  wall 
layer  and  3)  a blending  layer  of  thickness  intermediate  between  that  of  the  outer 
and  wall  layers.  In  this  respect,  the  theory  is  superficially  similar  to  the 
"tr iple-de^k"  structure  arising  in  the  laminar  interaction  theory  of  Brown  and 
Stewar tson  . 

Although  they  both  envolve  a three-layer  structure,  the  details  of  the 
laminar  and  turbulent  solutions  are  vastly  different.  In  the  laminar  problem 
the  interaction  distance  is  an  order  of  magnitude  larger  than  a boundary-layer 
thickness,  with  boundary-layer  approximations  and  the  equivalent  displacement- 
body  concent  remaining  valid  in  the  trailing -edge  region.  In  turbulent  flow,  the 
interactions  are  compact  with  a streamwise  extent  of  the  order  of  a boundary- 
layer  thickness,  the  boundary-layer  approximation  fails  and  the  displacement 
surface  is  not  a streamline  in  the  trailing-edge  region.  The  reasons  for  these 
differences  will  be  discussed  during  the  talk. 
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W'i!  will  also  discuss  how  the  turbulent-interaction  equations  for  the  outer 
solution  are  automatically  closed  to  lowest  order  without  the  need  for  empirical 
closure  assumptions.  We  will  outline  the  procedures  for  obtaining  solution  of 
the  inviscid  outer  equations  by  analytic  function  theory.  The  effect  of 
compressibility  on  the  solution  will  be  discussed  and  a^compiete  solution  for  the 
lift  coefficient  in  a turbulent,  compressible  flow  will  be  given. 
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FLOW  SEPARATION  FROM  A SMOOTH  SURFACE  IN  A SLENDER  CONICAL  FLOW 


N.  Riley 

University  of  East  Anglia,  Norwich,  U.K. 


ABSTRACT 


When  a flat  plate  delta  wing  is  placed  at  incidence  to  an  oncoming  stream 
the  dominant  feature  of  the  resulting  flow  is  the  vortex  pair  which  forms  over 
the  wing,  as  in  figure  1.  The  flow  separates  from  each  leading  edge,  Sj,  to 
feed  vorticity  into  these  coherent  vortices.  Efficient  calculation  methods 
are  avai lable 1 » --  for  the  investigation  of  this,  so-called,  primary  separation 
for  a slender  wing.  It  is  known  that  the  pressure  variation  over  the  wing 
surface  due  to  the  presence  of  these  leading-edge  vortices  induces  a secondary 
separation.  It  seems  natural  to  explore  the  possibility  of  representing  this 
secondary  separation,  like  the  primary  separated  flow,  by  a conical  spiral 
vortex.  The  model  is  shown  schematically  in  a cross-flow  plane  in  figure  2, 
where  S2  denotes  the  position  of  the  secondary  separation  line.  Numerical 
calculations  have  been  carried  out-^,  in  the  boundary  layer  which  forms  on  the 
surface,  to  determine  the  position  of  S2  such  that  boundary-layer  separation 
also  occurs  at  that  point.  It  has  been  found  that  separation  takes  place  at 
or  just  beyond  a pressure  minimum  in  this  conical  flow.  This  is  reminiscent 
of  recent  theories  for  two-dimensional  separation^.  A local  investigation  of 
the  inviscid  flow  structure  close  to  S2  has  been  made-’  which  shows  a close 
resemblance  to  the  discontinuous  potential  flow  models  of  separation  in  two 
dimensions.  The  main  difference  is  that  in  the  present  case  the  region  beyond 
St  is  not  a stagnant  region.  An  investigation  of  the  viscous  interaction  flow 
structure  is  currently  in  progress, exploi ting  the  ideas  which  have  been 
developed  for  two-dimensional  flow^*. 


O - O - 

1 Figure  1 1 Figure  2 
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DIFFICULTIES  WITH  BRUTE-FORCE  NUMERICAL  SOLUTIONS  OF 
NAVIER-STOKES  EQUATIONS  FOR  HIGH  - RE  SEPARATED  FLOWS 

Patrick  J.  Roache 

Albuquerque,  New  Mexico  87108 

ABSTRACT 


Discrete  numerical  solutions  in  fluid  dynamics  are  well-suited  for 
problems  in  an  intermediate  range  of  Re,  which  are  very  difficult  to  approach 
by  analytical  methods.  At  high  Re,  difficulties  are  encountered  which  were 
generally  unanticipated  when  computational  fluid  dynamics  was  in  its  infancy. 

These  difficulties  will  be  viewed  from  several  different  perspectives: 
"balancing"  a method;  the  "swamping"  of  viscous  terms  by  truncation  error  from 
advection  terms;  spatial  oscillations;  and  most  importantly,  scaling  and  the 
Nyquist  frequency  limitation.  These  perspectives  give  estimates  of  the  maximum 
Re  resolvable  in  a given  mesh.  However,  computational  experience  indicates  that 
these  estimates  are  pessimistic  if  only  moderate  accuracy  is  desired.  On  the 
other  hand,  success  with  a moderately  accurate  solution  can  be  misleading  about 
the  prospects  for  high  accuracy  solutions.  If  a problem  can  be  solved  with 
second-order  methods,  it  can  be  solved  to  higher  accuracy  with  high-order 
methods,  but  high-order  methods  do  not  help  the  scaling  problem  at  high  Re. 

Perspectives  will  also  be  given  on  what  constitutes  acceptable  accuracy 
for  some  applications,  with  examples  of  calculations  which  are  accurate  for 
some  parameters  (e.g.  drag)  but  not  for  others  (e.g.  pressure  distribution), 
or  which  give  an  "acceptable"  error  at  one  Re,  but  do  not  represent  variation 
with  Re. 

Some  possible  techniques  for  obtaining  acceptable  accuracy  will  be  reviewed, 
including  the  following:  (l)a  judicious  choice  of  problems,  which  sometimes 
means  that  the  problem  would  be  more  efficiently  calculated  with  some  form  of 
parabolized  Navier-Stokes  equations,  (2)  for  rough  accuracy,  the  use  of  implicit 
or  explicit  "artificial  viscosity",  which  is  directional  and  is  better  than 
just  calculating  at  a lower  Re,  (3)  automatic  mesh  refinement  and  coordinate 
alignment,  (4)  possibly,  the  inverse  approach  of  Harlow  for  one-dimensional 
problems,  which  is  conceptually  similar  to  the  isotherm  migration  technique 
used  in  heat  conduction  problems,  (5)  possibly,  the  incorporation  of  results 
from  perturbation  analyses. 
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THE  CALCULATION  METHODS  FOR  FLOWS  WITH  STRONG  INTERACTIONS 
BETWEEN  VISCOUS  AND  INVISCID  REGIONS 
J.  Rom 

Technion-  Israel  Institute  of  Technology 
Haifa,  Israel 

ABSTRACT 


In  flows  with  strong  interaction,  particularly  when  separated  regions 
are  present,  small  changes  in  the  viscous  layer,  and  conversely  small  dis- 
turbances in  the  external  pressure  distribution  may  precipitate  large  dis- 
tortions of  the  viscous  layer  and  in  the  flow  field.  In  such  cases  the  most 
generally  valid  approach  would  be  a numerical  solution  of  the  full  Navier-Stokes 
equations,  although  such  solutions  are  hampered  by  the  computer  size  and  speed 
limitations1.  Therefore,  it  is  recognized  that  there  exists  the  need  for  a 
theory  capable  of  giving  good  representation  of  the  flow  field  in  the  regions 
of  strong  interaction,  without  developing  singular  Icehaviour  at  such  points  as 
separation.  Such  flow  fields  are  very  difficult  to  conpute  because  it  is 
found  that  equations  of  hyperbolic,  parabolic  and  elliptic  types  enter  the 
problem  in  different  regions  and  have  to  be  joined  smoothly  on  the  interface 
between  neighbouring  regions.  Furthermore,  the  characteristics  of  the  viscous 
and  inviscid  regions  and  the  compatibility  relations  on  the  matching  line  result 
in  an  additional  scaling  parameter  for  the  viscous  flow'1.  An  experimental 
result  illustrating  this  effect  is  shown  in  Fig.  1.  This  means  that  in  scaling 
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Fig.  Is  Variation  of  separation  length  with  Reynolds  number (from  Jones  - 
AGARD  CP-83,  1971. 

experimental  results,  the  external  flow  Much  and  Reynolds  numborfmust  be  matched  as 
well  as  the  relative  dimension  of  the  viscous  region  (including  any  separated 
zones).  In  tile  separated  zone,  the  proper  flow  reversal  velocities  can  be  obtained 
by  a solution  of  the  Navier-Stokes  equations.  This  will  result  in  a three  region 
model:  (1)  external  flow  (inviscid  -hyperbolic,  elliptic  or  mixed) , ( J)  shear 

region  (boundary  layer-parabolic),  (3)  separated  region  (flow  reversal-elliptic). 
In  most  cases  there  will  be  an  overlap  between  regions  (2)  and  (3)  and  the  com- 
patibility relations  will  be  satisfied  in  this  overlap  zone. 

The  solution  of  the  complete  flow  field  will  be  obtained  by  an  iterative 
procedure.  The  key  problem  for  convergence  of  such  iterative  procedures  are  the 
proper  matching  methods.  The  matching  can  be  done  either  by  the  schematic  pro- 
cedure suggested  in  Ref.  3,  or  by  a "Bisection  procedure"  defined  in  the 
following:  (1)  Assume  an  initial  velocity  distribution  on  the  outer  edge  of  the 

viscous  layer  UB{Xp},  (2)  Solve  by  a finite  difference  program  the  boundary  layer 
equations  for  the  imposed  UB{Xi)and  obtain  6B{X;}.  (3)  Assume  SgfXi  j*=6p{  Xi } , 
calculate  the  external  velocity  distribution  Up(x^).  (4)  Define  a new  v<l  :ity  for 
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the  outer  edge  of  the  boundary  layer  UgN{Xi}=  — [Un{Xi}+Up{Xi} ] . (5)  Repeat 

iterative  calculation  of  steps  (1)  to  (4)  with  UB{Xi}=UBN{Xi} . This  procedure 
is  simpler  than  that  of  Ref.  3 in  that  it  does  not  require  the  evaluation  of 
the  influence  matrices  [Bj_  j ] and  [Pi  j ] - Furthermore,  in  cases  where  the  bi- 
section procedure  does  not  converge,  that  is  when, 

. 3+1  j+l  j j 

|up  {X.}  - Ufi  {x.}|>|up  {X.}  - uB{x.}|, 

j 


then  computation  procedures,  illustrated  in  the  following  sketches,  can  be 


These  procedures  are  now  being  applied  to  the  investigations  of  the  following 
problems:  (1)  strongly  interacting  flow  in  the  vicinity  of  a forward  facing 

step  in  supersonic  flow.  (2)  calculation  of  the  interaction  of  a turbulent 
boundary  layer  in  a transonic  flow  on  a symmetric  airfoil. 
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VISCOUS  INTERACTIONS  IN  THREE  DIMENSIONS 


S . G . Rub i n 


Polytechnic  Institute  of  New  York 
Farmingdale,  New  York  11735,  U.S.A. 


ABSTRACT 


High  Reynolds  number  viscous  interactions  in  three-dimensions 
can  result  in  a variety  of  complex  flow  phenomena.  As  in  two- 
dimensions,  strong  pressure  interactions  and/or  mainstream  sep- 
aration may  be  important;  many  of  the  presentations  at  this  work- 
shop will  be  concerned  with  free  interaction  boundary  layer 
modelling  of  these  effects. 


Three-dimensional  interactions  introduce  the  possibility  of 
secondary  flow  separation  and  associated  vortex  formation,  as 
well  as  boundary  layer  interference  phenomena  associated  with 
geometry,  displacement  induced  secondary  motions  and  interactions 
near  lines  of  attachment  or  separation.  It  is  these  three- 
dimensional  viscous  interactions  that  have  been  of  interest  to 
me  for  sometime  and  about  which  this  presentation  is  devoted.  It 
should  be  noted  that  the  procedures  described  here  can  also  play 
an  important  role  in  the  analysis  of  free  interactions.  Investi- 
gations along  these  lines  are  currently  in  progress. 


Spycjj-fic  geometries  that  have  been  considered  by  merged  layer 
theory^)  (single  viscous -inviscid  system  of  equations),  two-layer 
theory^  ^youpled  viscous-inviscid  system)  boundary  region 
theory  ' ; (boundary,.  j.ayer  within  a boundary  layer)  , and  viscous 

slender  body  theory  ~ (second-oy^egj  evaluat^n^o^  ^ycondary 
motions) , ^yclude  axial  corners' ''  , cones'  ' , ogive- 

cylinders'  , elliptic  cylinders'  and  wing-body  configurations. 
Rectangular  inlets  have  also  been  evaluated  but  only  external 
interactions  are  included  here.  Our  main  objective  in  all  of 
these  analyses  is  the  evaluation  of  complex  flow  phenomena 
without  requiring  complete  solutions  of  the  Navier-Stokes  equa- 
tions. Some  typical  secondary  motions  resulting  from  these 
analyses  are  given  below.  Figure  1 shows  the  separated  secondary 
flow  profile  (g)at  the  leeplane  of  a cone  at  incidence; in 
reference  7 the  inadequacy  of  conventional  boundary  layer  methods 
is  clearly  defined  and  the  odd  growth  for  large  £ is  explained. 
Figure  2 depicts  the  low  speed  secondary  flow  on  a cylinder 
resembling  a wind-body  cross-section.  This  is  a recent  result 
obtained  with  viscous  slender  body  methods. 
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Fig.  1,  Similar  Boundary  Layer:  Velocity  Profile  at  Leeplane. 


Fig.  2.  Secondary  Flow  on  Wing-Body  at  Small  Incidence. 
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Tube— flows  : on  the  effects  of  constrictions,  branchings, 
upstream  influence  and  separation 


F.T.  Smith 

Mathematics  Department,  Imperial  College,  London  S.W.7*  England 

ABSTRACT 


These  tubeflows  are  divided  into  four  main  categories:  (i)  non symmetric 

pipeflows  (ii)  nonsymmetric  channel  flows  (iii  ) symmetric  pipeflows 
(iv)  symmetric  channel  flows.  The  Reynolds  number  R is  large  and  the 
incompressible  fliud  motion  is  steady  and  laminar.  The  aim  is  to  predict 
flow  features  in  the  most  realistic  case  where  the  obstruction  (a 
constriction  or  branching  of  the  tube)  is  ''severe'',  i.e.  of  a size 
comparable  with  the  undisturbed  tubewidth  (a). 


Flow  Cii)  provides  the  most  upstream  influence  then,  over  a large  distance 
0(aHV7)t  and  upstream  separation  occurs  at  x = X-)Sp  = -0.49  a R1/?  + 0(a), 
independent ly  of  the  particular  form  of  obstruction.  (Here  x is 
distance  from  the  obstruction).  The  two  wall  layers  interact  with  the 
core,  one  layer  expanding,  while  the  other  compresses  and  s ■ sustains 
the  upstream  interaction  process.  Comparisons  with  experiments  and  with 
solutions  Of  the  Navien-Stokes  equations  are  favourable  overall. 


2 

Flows  (iii),  (iv)1-  are  treated  together  and  thi  upstream  infiueno 
li stance  (and  hence  — x„Cp)  is  0(a  lnR).  For,  in"r  * ni  • rictions 

of  slope  h R“l/6  (h  /vj  l1),  xsep  is  0(l)  but  as  h — xsep  • - oo  . 

Hence  in  severe  obstruction  x0  is  again  independent  f the  obstruction 
shape,  to  first  order. 

Flow  (i)’  is  the  most  physically  important ,but  as  yet  we  can  only  conclude 
•hat,  for  a ''fine'1  contriction  (of  sloped  R~ 1/3;,  upstream  influence 
exists  over  a distance  0(a).  Th  -operties  for  severe  obstructions 
.•wait  a numerical  investigation  of  the  complete  3D  free— interaction  that 
occurs  for  the  fine  case.  This  interaction  is  due  to  the  jetlike  swirl 
in  the  wall  layer,  in  the  absence  of  any  core  displacement. 
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SYCHEV'S  THEORY  OF  SEPARATION 

K.  Stewartson 

The  Ohio  State  University 
Columbus,  Ohio  43210 

and 

University  College,  London  WCIE  6BT,  U.K. 


ABSTRACT 


The  limit  of  the  solution  of  the  Navier-Stok.es  equations  for 
flow  past  a bluff  boay  as  the  Reynolds  number  R -*•  00  is  not  yet 
known.  Two  candidates  have  been  sufficiently  well  described  that 
one  may  ask  whether  they  are  the  leading  term,  not  uniformly  valid, 
in  an  asymptotic  expansion  of  the  solution  in  descending  powers  of  R. 
The  first,  lrrotational  attached  flow,  seems  to  fail  this  test.  The 
second  is  Kirchhoff  flow  of  which  the  most  likely  has  breakaway  at  a 
point  S of  the  surface  such  that  the  free  stream  line  there  has  zero 
curvature  and  the  pressure  gradient  upstream  is  favourable.  By  the 
use  of  a triple-deck  theory  to  describe  the  boundary  layer  near  S, 
Sychev  establishes  that  the  description  is  viable  here  provided  a 
certain  differential  equation,  of  the  boundary  layer  type,  has  a 
solution.  Some  comments  are  made  about  this  equation  and  the  relation 
of  the  theory  to  existing  studies  of  flow  past  bluff  bodies  at  finite 
values  of  R. 

The  properties  of  flow  past  finite  thin  bodies  at  incidence  when 
R >>  1 is  reviewed,  it  being  noted  that  complete  solutions  have  been 
found  up  to  but  not  beyond  separation.  Possibilities  for  solutions 
with  separation  are  touched  on. 

Sychev  went  on  to  claim  that  the  notion  of  separation  being 
preceded  by  a finite  length  of  adverse  pressure  gradient  has  no 
place  in  a theory  of  flow  past  a bluff  body  when  R >>  1.  The 
strength  of  the  case  for  the  notion  in  laminar  flow  is  reviewed  and 
the  implications  of  his  claim,  if  substantiated,  is  discussed  for 
both  laminar  and  turbulent  flows. 
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VISCOUS  FLOW  NEAR  THE  TRAILING  EDGE  OF  A FLAT  PLATE 


A. I.  van  de  Vooren 

Mathematical  Inst  itute.  Uni  versity  of  Groningen 
Netherlands 

Re  earch,  performed  during  the  last  year  at  the  Mathematical  Institute  of  the 
University  of  Groningen,  mainly  by  dr.  Veldmarv  will  be  reviewed. 

The  equations  of  the  triple  deck,  existing  near  the  trailing  edge  of  a flat 
plate  under  zero  angle  of  incidence,  have  been  solved  numerically.  By  aid  of 
suitable  transformations  the  whole  range  of  X-values  (from  to  +*»  on  the 
triple  deck  scale  X=ReJ'°x)  has  been  considered  without  cutting-off.  Another 
difference  with  other  investigator's  work  is  that  the  solution  in  the  upper 
deck  has  been  obtained  by  the  aid  of  a field  calculation  which  is  performed 
simultaneously  with  the  field  calculation  of  the  lower  deck.  This  prevents  the 
use  of  the  Cauchy  integral  whcih  either  leads  to  difficulties  at  X=0  or  at 
X=±  <*>.  In  the  usual  boundary  layer  and  in  the  wake  (|x|  OfRe-^^))  the  displace- 
ment thickness  determines  the  pressure  distribution,  while  in  the  regions  inside 
the  triple  deck  ( j xj  - 0(Re~3/8) ) it  is  the  (self-induced)  pressure  which  determines 
the  displacement  thickness.  Therefore,  for  x=0(Re_i,/®)  a transition  from  the  one 
to  the  other  system  occurs  and  this  makes  that  the  direction,  in  which  the 
Cauchy  integral  should  be  applied,  also  changes. 

On  an  x-scale  smaller  than  that  of  the  usual  triple  deck  an  infinite  series 
of  triple  decks  occur.  Their  sizes  are  x=0(Re-' /d) , x=0(Re-^/®) , etc.,  finally 
converging  to  the  Navier-Stokes  region,  x=0(Re-3/4).  However,  the  main  term  in 
the  shear  stress  at  the  wall  remains  constant  from  the  inner  end  of  the  orthodox 
triple  deck  until  the  outer  end  of  the  Navier-Stokes  region.  The  discontinuity 
in  the  pressure  derivative  at  X=0  disappears  in  the  upper  layer  of  the  new 
triple  decks.  It  remains  present  in  the  middle  and  lower  layers  of  these  triple 
decks  until  the  Navier-Stokes  region  is  reached. 

For  X-'  the  expansions  of  stream  function  und  pressure  in  the  wake  contain  so- 
called  eigenfunctions.  Their  coefficients  can  not  be  determined  from  the  asymp- 
totic theory  but  only  follow  from  the  numerical  determination  of  the  whole  flow 
field. 

Finally,  a momentum  consideration  is  given  which  leads  to  a relation 

d=-6  a b 

where  d is  the  coefficient  of  the  Re'^/^-term  of  the  drag,  a = f''(0)  the  second 
derivative  of  the  Blasius  profile  at  the  wall  and  b is  the  coefficient  of  the 
most  important  eigenfunction. 


NUMERICAL  STUDIES  OF  HIGH  REYNOLDS  NUMBER  LAMINAR  SEPARATION 

M.J.  Werle  and  R.T.  Davis 

University  of  Cincinnati 
Cincinnati,  Ohio,  45221,  U.S.A. 

ABSTRACT 

The  goal  of  our  ongoing  effort  is  to  help  establish  a fundamen- 
tal set  of  governing  equations  and  a solution  technique  for  reliable 
prediction  of  separation.  The  approach  taken  attempts  to  comprehend 
as  much  analytical  information  as  possible  in  the  numerical  solution 
method.  Results  are  discussed  for  supersonic  ramp  and  injection 
induced  separations  and  for  incompressible  flow  separation  at  ex- 
terior (convex)  and  interior  (concave)  corners. 

For  supersonic  induced  separation,  Fig.  1 shows  that  numerical 
solutions  of  the  triple  deck  equations  (TD)  predict  the  overall 
scale  of  separation  as  obtained  from  numerical  solutions  of  the  in- 
teracting boundary  layer  (IBL)  equations  as  Re-*°°.  Most  importantly 
it  was  found  that  accurate  solutions  to  the  IBL  equations  could  not 
be  achieved  unless  the  finite  difference  mesh  honored  the  scalings 
predicted  by  the  TD  model.  Also  it  has  been  found  that  the  parame- 
ter groupings  that  emerge  from  the  supersonic/hypersonic  TD  model 
correlate  IBL  separation  solutions  thereby  reducing  the  4 fold  de- 
pendence on  M^,  Tw,  Re^  and  a to  a 2 parameter  family  - a result 
that  should  be  helpful  to  design  studies.  Fig.  1 also  shows  results 
of  an  integral  method  (IM)  solution  of  the  TD  equations  obtained 
using  a quadratic  velocity  profile.  The  simplicity  of  the  IM  ap- 
proach makes  it  attractive  not  only  as  a simple  solution  method  for 
studying  the  fundamental  interaction  process  but  also  as  a vehicle 
for  obtaining  (a)  initial  conditions  for  relaxation  solutions  of 
the  IBL  model,  (b)  the  overall  size  of  the  interaction/separation 
flow  region,  and  (c)  a local  slip-boundary  condition  for  more 
elaborate  solution  schemes.  For  the  ramp  problem  there  are  several 
areas  for  future  study  - these  include  development  of  a means  for 
accommodating  imbedded  shock  waves,  a better  understanding  of  the 
mechanism  for  upstream  propagation  of  information,  and  a faster  re- 
laxation scheme  for  solving  the  IBL  equations. 

Figure  2 shows  that  application  of  the  IBL  model  to  slot  injec- 
tion into  a supersonic  stream  verifies  the  TD  prediction  of  separa- 
tion occurring  ahead  of  the  slot  and  verifies  that  "blow  off"  will 
not  occur.  IBL  plate  injection  studies  indicate  the  same  trend  but 
such  solutions  were  difficult  to  obtain  because  of  the  build  up  of 
an  inviscid  layer  near  the  wall  and  a slow  convergence  to  an  asymp- 
totic downstream  flow  state.  Further  study  is  needed  of  the  numer- 
ical method  for  such  viscid/inviscid  flows  and  ultimately  correla- 
tions based  on  TD  concepts  should  be  developed. 

For  incompressible  flow  separations,  numerical  studies  have 
been  conducted  of  the  flow  past  interior  and  exterior  corners  on 
semi-infinite  bodies  by  mapping  the  geometries  to  stagnation  point 
flows.  A vorticity  stream  function  formulation  of  the  IBL  concept 
has  been  followed  with  typical  wall  shear  distributions  for  the  ex- 
terior corner  flow  shown  in  Fig.  3.  It  is  found  that  for  low  Re 
solutions  are  easily  obtained  whereas  for  high  Re  the  solutions  are 
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Streaming  Motion  Past  Bluff  Bodies  at  Re  from  0(1)  to  0(10?) 
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ABSTRACT 

A new  approximate  theory  will  be  described  for  treating  the  flow  past 
smoothly  contoured  two-dimensional^  and  axisymmetric^  bluff  bodies  in  the 
intermedia^0  Reynolds  number  range  0(1)  < Re  < 0(1 02 ) where  the  displacement 
effect  of  the  thick  viscous  layer  near  the  surface  of  the  body  is  large  and  a 
steady  laminar  wake  is  present.  The  theory  is  based  on  a new  pressure  hypo- 
thesis which  enables  one  to  take  account  of  the  displacement  interaction  and 
centrifugal  effects  in  thick  viscous  layers  using  conventional  first  order 
boundary  layer  equations.  The  basic  question  asked  is  how  the  wall  pressure 
gradient  in  ordinary  boundary  layer  theory  must  be  modified  if  the  pressure 
gradient  along  the  displacement  surface  using  the  Prandtl  pressure  hypothesis 
is  to  be  equal  to  the  pressure  gradient  along  this  surface  using  a higher  order 
approximation  to  the  Na'M'er-Stokes  equation  in  which  centrifugal  forces  are  con- 
sidered. The  inclusion  of  the  normal  pressure  field  with  displacement  inter- 
action is  shown  to  be  equivalent  to  stretching  the  streamwise  body  coordinate  in 
first  order  boundary  layer  theory.  For  laminar  flow  without  separation  the  new 
theory  provides  a smooth  transition  to  the  solutions  obtained  from  high  Reynolds 
number  boundary  layer  theory. 

While  the  new  theory  is  of  a non-rigorous  nature,  it  yields  results  for  the 
location  of  separation  and  detailed  surface  pressure  and  vorticity  distribution 
which  are  in  remarkably  good  agreement  with  the  large  body  of  available  numeri- 
cal Navier-Stokes  solutions.  These  solutions  include  the  flow  past  parabolic, 
circular  and  elliptic  cylinders,  spheres,  speroids  and  paraboloids  of  revolution 
at  various  Reynolds  numbers.  A novel  feature  of  the  new  boundary  value  problem 
is  the  development  of  a simple  but  accurate  approximate  method  for  determining 
the  inviscid  flow  past  an  arbitrary  two-dimensional  or  axi symmetric  displacement 
body  with  its  wake. 
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A METHOD  FOR  REVERSE  FLOW  BOUNDARY  LAYERS 
WITH  PRESCRIBED  DISPLACEMENT  THICKNESS 
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ABSTRACT 


a method  is  described  for  computing  boundary  layers  with  pre- 
scribed displacement  thickness,  together  with  the  modifications 
needed  to  implement  the  iterative  procedure  introduced  in  Ref.(l) 
for  dealing  with  reverse  flow  regions.  The  main  aim  is  to  inves- 
tigate more  fully  the  convergence  of  this  iterative  procedure  on  a 
simpler  proolem,  namely  an  incompressible  boundary  layer  with  dis- 
placement thickness  prescribed  in  such  a way  that  a bounded  reverse 
flow  region  is  obtained.  This  avoids  the  problem  of  having  to  im- 
pose a downstream  asymptotic  profile,  which  was  necessary  in  Ref.(l) 
where  the  convergence  of  the  iterative  procedure  was  not  completely 
convincing.  It  should  be  remarked  that  this  may  have  been  due  to 
a low  order  term  in  the  asymptotic  expansion,  given  in  Ref. (2), 
being  incorrect,  as  has  been  kindly  pointed  out  to  us  by  Dr.  F.T. 
Smith;  the  effect  snould  be  small,  but  whether  it  actually  removes 
the  slight  difficulty  noted  in  Ref.(l)  has  not  yet  been  confirmed. 


The  specific  problem  considered  here  has  a displacement  thick- 
ness given  by 


a*(x) 


<5  0 ( 1 + 9e 


-(x-14)2/16} 


wnere  6-  is  the  constant  displacement  thickness  of  the  stagnation 
point  boundary  layer.  Thus  we  have  basically  a stagnation  point 
boundary  layer  modified  in  such  a way  that  the  displacement  thick- 
ness nas  a hump  centered  at  x=l4  but  drops  smoothly  to  the  constant 
6q  for  large  x and  also  for  small  x for  practical  purposes.  Actu- 
ally it  was  assumed  that  a stagnation  point  boundary  layer  had 
developed  up  to  x=l  and  was  then  allowed  to  develop  in  accordance 
with  prescription  of  6*(x). 

The  basic  step-by-step  scheme  used  is  that  of  Ref.(l)  with 
modifications  to  account  for  the  different  boundary  conditions.  In 
Ref.(l)  the  inner  layer  of  a compress! b le  interaction  problem  was 
treated  and  the  outer  boundary  conditions  were  u-*y-A(x),  where 
A=  p dx,  and  Uy-*-l.  These  are  here  replaced  by  u-*-U(x),  where 
P+  1/2  U2  = constant,  and  /q  ( U-u ) dv*6 * ( x ) . We  may  remark  that  the 
method  could  also  be  modified  for  the  corresponding  inner  layer 
problem,  where  A(x)  is  specified  Instead  of  the  relation  between 
A and  p,  which  may  be  used  in  treating  incompressible  interaction 
problems  in  an  iterative  manner. 

The  procedure  for  the  reverse  flow  region  begins  by  producing 
a first  approximation  by  the  technique  first  described  by  Fliigge-Lot 


and  heyhner3  as  referred  in  Hef.(l)  as  FLARE.  This  is  then  corrected 
successively  by  downstream-upstream  iterations  referred  to  as  DUITS. 
Eacn  upstream  iteration  starts  at  the  end  of  the  reverse  flow  bub- 
ble and  is  confined  to  the  reverse  flow  region  produced  by  the 
previous  downstream  iteration  and  assumes  the  pressure  distribution 
obtained  then  as  given.  Each  downstream  iteration  determines  the 
pressure  to  satisfy  the  given  <5#(x)  and  in  the  reverse  flow  region 
takes  the  uux  term  as  given  by  the  previous  upstream  iteration; 
the  first  approximation  is  obtained  by  setting  this  equal  to  zero. 

A maximum  of  eight  DUITS  were  performed  and  convergence  appeared 
to  be  satisfactory.  A reasonable  practical  accuracy  was  obtained 
after  three  or  four  DUITS. 
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